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Synopsis. 


The problem of the alpha fine structure intensities for deformed nuclei is 


treated on the basis of the unified nuclear model. A system of coupled diflerential 
equations for the radial wave functions is obtained; they correspond to the dif- 
ferent channels into which the alpha particle is able to leak out, leaving the daughter 
nucleus in different final states. An approximate analytical formula is derived for 
the intensities of the different alpha groups in terms of the alpha wave function 
on the nuclear surface. Also an analysis of the empirical data is presented. 

The nuclear deformation causes the alpha decay rate to be enhanced by a 
factor which increases slowly with the nuclear deformation and which may amount 
to about 2 for the largest deformations actually found. However, an enhancement 
factor of this kind is difficult to detect experimentally. 

The empirical intensities of the even parity alpha groups of even-even nuclei 
indicate the existence of higher multipole deformations of even parity in the 
nuclear surface. The magnitudes of these higher order deformations may change 
rapidly with the atomic number Z. From the empirical intensities of the odd parity 
alpha groups which have been found in some even-even nuclei, an estimate is 
presented of the magnitude of the nuclear octupole moment, the existence of 
which has been suggested by CHRrisryY to explain the presence of the odd parity 
alpha groups for these even-even nuclei. 

The favoured alpha decays of odd-A nuclei and even-even nuclei are similar, 
and it is possible to obtain an approximate formula which relates the intensities 
of the favoured alpha groups of an odd-A nucleus to the intensities of the alpha 
groups of neighbouring even-even nuclei. This intensity formula which was first 
given in a recent paper by Bonr, FROMAN, and MOTTELSON is discussed in greater 
detail. 

The angular distributions of alpha particles from polarized nuclei are also 
considered. 


PRINTED IN DENMARK 
BIANCO LUNOS BOGTRYKKERT A/S 


I. Introduction. 


few years after the advent of quantum mechanics, GAMOW, CONDON, and GURNEY 

explained the puzzling features of alpha decay then known by treating the alpha 
decay as a barrier penetration problem. (See, e. g. GaMow and CRITCHFIELD?8,) In 
this way, they were able to account approximately for the dependence of the half-life 
on the alpha particle energy and on the atomic number of the alpha emitting nucleus. 
Shortly afterwards, ROSENBLUM discovered the alpha fine structure which is due to 
the fact that the alpha particle may leave the daughter nucleus in different low-lying 
states. At that time, one could not account for the details of the fine structure pattern 
and for the intensities of the alpha groups. 

During the last decade, extensive and accurate experimental work has been done 
in alpha spectroscopy, which has revealed many striking regularities in the alpha 
fine structure pattern. (See the review article by PERLMAN and Asaro??,) These re- 
gularities have been found to be related to the rotational level structure which is 
characteristic of nuclei possessing a shape differing strongly from spherical symmetry. 
Such large nuclear deformations are known to occur for nuclei with configurations 
far removed from closed shells, as, for example, in the region of heavy elements 
(4 >220). (See BoxHr and MotTELson?"f,) 

For the interpretation of the alpha fine structure, it is therefore imperative to 
take into account the departure of the nuclear shape from spherical symmetry. The 
nuclear deformation is important not only for the energies of the different alpha 
groups, but also for their relative intensities. Furthermore, as was first pointed out 
by Hixrx and WHEELER?, the nuclear deformation may also affect both the absolute 
alpha decay rate and the angular distribution of the alpha particles. 

In a previous paper by Boxr, FRromAN, and MotTELsonN'””, in which the inter- 
pretation of the alpha decay systematics is discussed on the basis of the unified nuclear 
model, certain simple relations between the relative intensities of the alpha groups 
are found. In the present paper, the problem of the alpha decay of deformed nuclei 
is treated in more detail and this treatment is also coordinated with a new analysis 
of the experimental data. In particular, we derive approximate formulae for the in- 
tensities of the different alpha groups in terms of the alpha wave function on the 
nuclear surface and the shape of this surface. It is found that the intensities depend 
very sensitively not only on the quadrupole moment, but also on higher multipole 
moments of the nuclear shape. Thus, the observed alpha intensities may yield in- 
formation on the details of the nuclear shape as well as on the probability for alpha 
particle formation in different parts of the nuclear surface. 

1= 
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In Chapter II, we give a brief review of the description of strongly deformed 
nuclei according to the unified nuclear model. The treatment of the alpha penetration 
problem is discussed in Chapters III—VI. Chapter VII contains an analysis of the 
experimental data. The applicability of the theory to the account of empirical data 
is discussed in Chapter VIII. The paper is followed by two appendices (and a list 
of references is added). Appendix A contains formulae and curves for the Gamow 
wave functions which describe the penetration of alpha particles through a spherical 
Coulomb barrier. Appendix B deals with the penetration of an alpha particle through 
an anisotropic potential barrier which is fixed in space. This problem is treated by 
combining a three-dimensional WKB-method, first suggested by CHristy, with ideas 
known from optics. 

This work was begun in 1953—1954 when the author was a member of the 
Theoretical Study Group of CERN (Conseil Européen pour la Recherche Nucléaire) 
at Universitetets Institut for Teoretisk Fysik in Copenhagen. I wish to thank Professor 
NIeLs BOHR for the privilege of working at his institute. The present investigation 
was suggested to me by Professor AAGE BoHR, to whom I am deeply indebted for 
many valuable discussions and advice, and for his continued interest in this work. 
I am grateful to Dr. BEN MOTTELSON for fruitful discussions of the problem and to Pro- 
fessor I. PERLMAN from the University of California, Berkeley, for numerous data put 
at my disposal previous to publication. To Professor Ivar WALLER, University of Upp- 
sala, I wish to express my gratitude for what he has taught me during the time when 
I was his student as well as for his encouragement in the course of this work. Finally, 
it is a pleasure to thank Mrs. S. HELLMANN for revising the English text and for her 
untiring help with the correction of the proofs. 


II. Survey of Some Features of the Unified Nuclear Model. 


The systematic occurrence of rotational nuclear spectra in the region of the very 
heavy elements shows that it is possible for these nuclei to separate approximately 
between an intrinsic motion of the nucleons in a deformed field of axial symmetry 
and a collective rotational motion of the nucleons. Corresponding to this separability 
of the motion, the normalized nuclear wave functions are; 


omen ' Sø | 
Wu, Sk re OD (IT-1) 
11k 0 and 
Ym,0—V/21+1 %xDi, 9 (&) (1-2) 


if K=0. In these expressions, the functions %,,, which are assumed to form an 
orthonormal set of functions describe the intrinsic particle SBÅ vibrational motion 


£ The normalization factors in (II-1) and (II-2) differ by a factor (8 Pe) from the corresponding 
normalization factors used by BoHR and MorrteELson-f, owing to the fact that we define the volume element 
dw such that id =1, whereas BoHRk and MOTTELSON use the definition | dw = SE 
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characterized by the quantum numbers specifying the particle configuration and the 
excited phonons in the nuclear vibrations. The quantum number K represents the 
component of the total angular momentum along the nuclear symmetry axis. If nothing 
else is stated, K is assumed to be non-negative. It seems possible to consider the in- 
trinsic particle motion approximately in terms of the independent motion of the 
individual nucleons, i. e., to consider the individual particle part of xx to be described 
by a Slater determinant depending on the space and spin coordinates of the individual 
nucleons relative to a coordinate system which is fixed with respect to the nuclear 
ellipsoid”, For the lowest rotational band of even-even nuclei, the nucleons move 
in paired orbits, and K = 0. In this case, the wave functions are given by (29): 
For odd-Å nuclei, odd-odd nuclei, and even-even nuclei with unpaired configurations 
K+ 0 and the wave functions are given by (II-1). Vibrational oscillations of the 
nuclear shape may be superimposed on the individual particle motion, but the ex- 
pected large excitation energy for such vibrations (Z 1 MeV) makes them relatively 
unimportant for the alpha decay problem. An exception may be the odd parity states 
which appear in some even-even nuclei and which may represent an excitation of 
a very soft asymmetric vibration. The symbol q in (II-1) is the sum of two quantities, 
one depending only on the intrinsic particle state, the other being the number of 
excited phonons corresponding to spherical harmonics of odd order in the nuclear 
vibrations. To specify clearly the nuclear vibrational state it is sometimes convenient 
fomusetthesgquantitytgkaskasfurtherfindexfoney ti est tor write ry instead dk OR Er 

The rotational nuclear motion is described by the D-functions which depend 
on the Eulerian angles w defining the orientation of the nucleus. These D-functions 
are characterized by the quantum numbers I, K, and M which are the total angular 
momentum of the nucleus, its component along the nuclear axis, and its component 
along a direction fixed in space (the z-axis), respectively. If the volume element dw 
is defined such that | dw=1, the system of functions V/21+1 DØ +Kk (w) is ortho- 
normal. The same is then also true for the system of wave functions p%, gf defined 
by (II-1) or (II-2). 

The energy levels of the nucleus are-given as the sum of one term, Tyart (K)", 
which depends on K and on other quantum numbers specifying the intrinsic state, 
and another term, TYy0t (K, I), which depends on both K and I. For a nucleus described 
by the wave functions (II-1), the possible values of I are 


LEE KE TE KER STES 


and the corresponding states of the nucleus have the same parity as the intrinsic state. 
Furthermore, Trot (K, I) in this case is given by 


te | mæ TALER | 
TER RED) = FN) NE NENS IE RE SEEGER SVED (Eg) 

«£ We note that this term may contain a contribution from vibrational as well as individual par- 
ticle motion. 
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where I is the effective moment of inertia, and ais the so-called decoupling parameter. 
For the ground state band of an even-even nucleus described by the wave functions 
(II-2), Tvart(0)=0 and the formula 


h2 
Trot (0, I) = 39 I(I+1) (II-4) 


applies. The possible values of I are 


KF 0332 AES 


and the corresponding states of the nucleus have even parity. Some of the very heavy 
even-even nuclei appear to have also a low-lying vibrational mode of excitation with 
K = 0 and odd parity. The nuclear wave functions are again given by (II-2), Tpart(0) 
has a magnitude of a few hundred keV, and TY%04(0, I) is given by (II-4) with I dif- 
fering somewhat from its value for the ground state band, and with the spin values 


TE TRES ekte 


The corresponding states of the nucleus have odd parity. 


III. Construction of the Wave Function of the Alpha Penetration Problem. 


The quantum numbers of the parent nucleus are recognized by an index i (for 
initial) and those of the daughter nucleus by an index f (for final). To fulfill auto- 
matically the conservation laws for the total angular momentum and for the angular 
momentum component along the z-axis, we describe the penetration of an alpha 
particle through the potential barrier by means of the wave function 


I; mt Er VENS pl EN 
UKE SER R > D le ETT) Dog (GEL Tp ON 0): (0) 
Se RE Er; 


where E is the total energy of the alpha particle and the daughter nucleus in the 
center-of-mass system, r is the distance between the alpha particle and the center 
of the daughter nucleus, and (9, w) are the polar angles of the alpha particle with 
respect to a center-of-mass coordinate system having axes with directions fixed in 
space. The Ø-functions appearing in (III-1) are defined by; 


db ; T 7 
Dy, (Ky,1, 14; w, 9, 9) = NE 13 Mein | 19,157 Mi) bx Yr my CO 9) SEE (15 
Å 
For the sake of simplicity, the dependence on the variables describing the intrinsic 
nuclear state has not been explicitly stated in this notation. lt is easily seen that the 
PD-functions form an orthonormal system of functions. The radial functions fx, BIG (T) 


£ All spherical harmonics and Clebsch-Gordan coefficients in this paper are defined according to 
CoNDON and SHORTLEY29, 
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in (III-1) are so far left undetermined. They depend on Kp, which is here used as 
an abbreviation for all the quantum numbers specifying the intrinsic state, on / and I, 
as well as on K; and I;. However, in the notation we have omitted, for the sake of 
simplicity, explicit reference to the dependence on K; and I;. The sum over Kyp which 
should be understood to include also a sum over all the other quantum numbers 
specifying the individual particle motion and the vibrational motion for the intrinsic 
state is included on the right-hand side of (III-1) in order to allow the description of 
transitions to different rotational bands of the daughter nucleus. 
The Clebsch-Gordan coefficients in (III-2) are different from zero only if 


Færn <=! < 1;+ I; 


and, therefore, the functions fra) are defined only if this inequality is fulfilled. 
This fact is also immediately evident from the rules for the addition of two angular 
momenta. The requirement of conservation of parity introduces the further restriction 
that / takes even values if the daughter nucleus is left in a band with the same parity 
as the parent nucleus, and that / takes odd values if the daughter nucleus is left in a 
band with a parity opposite to that of the parent nucleus. 

So far we have referred the motion of the alpha particle to a center-of-mass 
coordinate system with axes having directions, fixed in space, but this motion may 
also be referred to a center-of-mass coordinate system with axes having directions 
fixed with respect to the nuclear ellipsoid. AÅssuming this latter coordinate system 
to be a polar one with the polar axis along the axis of the nucleus, and denoting the 
coordinates of the alpha particle in this coordinate system by (r, &, p'), one realizes 
that it must be possible to transform the wave function (III-1) into the form 


De goes ED På FE 
mg 


IDI EN I; R RE, 
xV/: TE Fare Y,, QW, 2) 


iæ = IG y ' ' 
(31) 48% DUE e (0) Y,, QW P) 


(III-3) 


if "Kr +0, and into the form 


Ds RÅ y BEDS DE OG) VIEDE (0) Sr (ØS) (III4) 


21; FAR 


if Ky = 0. The functions kS= y in (III-3) and | IT; BF YND Mm, Q(&) VASER DN) 


in (III-4) form Mika søde systems of functions. 
The requirement that the appropriate one of the two wave functions (III-3) and 


(III-4) shall describe the alpha decay of a parent nucleus having the angular moment- 
DE: 
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um component K; along its axis shows that, on (or just outside) the nuclear surface, 
this wave function can contain only angular momentum components 2 fulfilling, 
KROER NEN KIG When is but slightly greater than the nuclear 
radius, the functions 9K4,1,2 (T) can therefore be different from 0 only for 2 = + K;—Ky. 
For K;=0, there exists only one such function if /> Ky, but none it EKK Or 
K;+0, there exist two if 1>K;+Kyp, one if |K:- K/|<I< K;+K/, but none if 
[= |K:-Ky|. Just outside the nuclear surface, the functions SKAL OD (7) are therefore 
different from zero for at most two Q-values, but, for larger values of r, these functions 
may be different from zero for all possible Q-values, since in general 2 is not a con- 
stant of the motion during the penetration (except for the decay of even-even nuclei). 
The derivation of the connection between the functions fx 1, TE) and IK, 1, BO) 
which we find in (III-9), requires some calculations. These calculations as well as 
those in the following chapter will be done for the case K/+ 0, i. e., assuming that 
the wave functions of the daughter nucleus are given by (II-1). However, in reality 
this is no restriction, for from (II-1) and (II-2) it is clear that the formulae for the 
case Ky = 0 can be obtained from the formulae for the case K+ 0 by leaving out 
the factor 1/y2 and the second term in the brackets i y, and by putting K/=0. 
By substituting (II-1) into (III-2), one finds 


I; 
Dj, (Kg,1,14 0,9, 9) | 
1 i. ren Es (III-5) 
= me UT; 0,20) HT (EN)? 1 X-k, PM (Kr 1,174; 0, &, »| | 
DÅ being defined by”; 
1 SER 
Er en 
j aA FENG (III-6) 
ST) VA 
= BE DI Ur, LÆNE m|Iy, kære Mi) V21;4+ 1 Dit, K,(W) GE CA P). | 


My Mm 
Using the formulae” 


Y,, m (&, 9) = » Dre), 0) 


Er l 
kr. TY NGE i É k 
SS (lg) GEM$: m | TAGER MER SSR SO TAI STEREO] Dig KRO.) 


k 
and 


VETTS KE OVE KEDE ED NOR ETGER Kr 40;=0 KLIT SKAN 
one easily finds that (III-6) becomesø 


£ In (III-6) and (III-7), as well as in the formulae between (III-6) and (III-7), the quantum number 
Kf is allowed to take also negative values. 


INT 3 9 


VEN CER gs 
Varg (RT, 00) =D] ERE TEE or) ve TE er | 
AED 

i I; ' ' 
xV27:+1D4,,x,+0 (2) Y1, 2 GF, 2). | 


From (III-5), (III-7), and the formula 
RER EU ET ER SE ENE ES ET EQT OVIDS RR) 
we then find that 


URET EL ØD 


I; 7 
KE 0,97) =D (21) 


CAD ORE O HIATT AE) 
ære 
AD 2 SE il TBA ' lå 
DE 50 2) (111-8) 
TI; 


TERE ED0 r r 
BEDE FE SEE SE CD) sgl 9) É 


This formula and the corresponding formula for the case Ky = 0 show that the wave 
function (III-1) can always be written in one of the forms (III-3) or (III-4) and that 
the connection between the functions (er re (7) and Te (7) is 


RTE Ø, ra : > 
dr. 00) == PD i Or ØE SØ | DEA CGETÆRKG ) fr enreC) (III-9) 
Si 


which can also be written as follows: 


. FREDE D Eg 7 
IraLT- (T) = DAG IEL ig (øsls Kr+Q,—Q|14, 1; 14, Ky) RAN OG) (III-10) 
For the alpha groups of even-even nuclei (1; = Ky = 0), these two formulae become 


VRD) (== 0) 00, 2 (III-11) 


fo rr.) ED) 6, 1p' (III-12) 


IV. Partial Transition Probabilities and Angular Distributions 
of the Alpha Groups. 


By an alpha group we denote those alpha particles which correspond to a tran- 
sition from a certain state (K;, I;) of the parent nucleus to a certain state (Kp, I4) of 
the daughter nucleus. It is the purpose of this chapter to give the formulae relating 
the partial transition probabilities and the angular distributions of the alpha groups 
to the asymptotic values (for large values of r) of the functions fk£,1,1, (7) which 
appear in the wave function discussed in the previous chapter. It is assumed that 
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the functions fky,1,14 (7) are normalized such that they correspond to the alpha decay 
of a single atomic nucleus. Furthermore, we consider the total energy EÉ as real, thus 
avoiding the unessential complication with complex eigenvalues and wave functions 
tending to infinity far away from the nucleus. 

Erom the wave function (III-1), which is valid both for K;+ 0 and for K/=0, 
one immediately finds that the probability per unit time, that an alpha particle leaks 
out with the angular momentum / leaving the daughter nucleus in a state (Ky, 14), is 


2. av 


BRA: = lim v|fx,,,1, 7) 


I taking either only even values or only odd values depending on the parities of the 
parent and daughter nuclei. The velocity of the alpha particle relative to the daughter 
nucleus is here denoted by v. The sum of the kinetic energy of the alpha particle 
and the kinetic recoil energy of the daughter nucleus is called the transition energy, 
and is denoted by Eftrans. The quantities Etrans and v are connected by the formula 


1 
Etrans = sm V23 


where m is the reduced mass of the alpha particle and the daughter nucleus. 
The partial half-life Tx4,1,7, which corresponds to the partial transition pro- 
bability Pr, RT, is 
In 2 


res (BIRDS 4 
LEGER ruste 


From (III-1) and (III-2) one easily finds that the probability per unit time and 
unit solid angle (in the center-of-mass system), that an alpha particle leaks out in 
the direction (%, p) leaving the daughter nucleus in the state (Kp, I, Mp), is 


2 
lim zu DICE 1; My, Mi — Mp 14,1 1, Må) fg ,4, 1407) 91, mom, Ø) (VS) 
r— 0 , 

where Or7,m (9) denotes the normalized Ø-depending factor of the spherical harmonic 
Yi,m (&, p). (See the book on the Theory of Atomic Spectra by CoNnDonN and SHORT- 
LEY20,) By summing the expression (IV-2) over the possible My-values and averaging 
it with respect to the initial distribution of M;-values, we can calculate the angular 
distributions of the different alpha groups emerging from polarized parent nuclei. 


V. Differential Equations of the Alpha Penetration Problem. 


In the center-of-mass system of the alpha particle and the daughter nucleus, 
the Hamiltonian of this two-body system is 


"ik 
H == ln Å,+ V (7) + Trart + fk. (V-1) 
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The reduced mass of the alpha particle and the daughter nucleus is denoted by m, 
the Laplacian with respect to r by 4,, the interaction energy between the alpha 
particle and the daughter nucleus by V (7), and the operator for the excitation energy 
of the daughter nucleus by Tpart+ Trot. The vectors r and r” are position vectors for 
the alpha particle relative to two different coordinate systems. The vector r = (r, d, %) 
is referred to a center-of-mass coordinate system whose axes have fixed directions 
in space, whereas the vector 7” = (r, %", p”) is referred to a center-of-mass coordinate 
system whose axes have fixed directions with respect to the nuclear ellipsoid of the 
daughter nucleus. 


The wave function BR which describes the hbarrier penetration of an alpha 


particle shall satisfy the wave equation 


h 0 il i 
(H+ EG) Bag 10 (V-2) 
and the boundary conditions that the alpha wave function is assumed to be known 
on the nuclear surface and that it shall represent only outgoing waves for large 
values of r. 

Starting from (V-1) and (V-2) we derive, in this chapter, differential equations 
for the functions fx, 11 (T) and IKy,1, 2 (Tr), which we have introduced in Chapter III. 
Although the calculations are here based on the wave function (III-3), i. e., we assume 
that K/+ 0, it is easily seen that the resulting differential equations for the functions 
LASE (T) and SA O (7) are valid also for K/;=0. 


For the Laplacian A, we have the well-known formula 


FE: 


me re 
7072 5 


where L is the orbital angular momentum (in units of hi). Expressed as a differential 
operator in the polar angles, L2 has, of course, the same form whether it is expressed 
in the angles (&, w) or in the angles (9, p"). From (III-3) and the above formula for 


ÅA, we get sliree 
TREE 8 ES TØS fl 
ÅR == (& h 222; dr r? Dit OLD) 
21;+1 I É reg . 
xJ/ SE | K/ PM, K/+ 2 (W) Yo (29) (VER) 


| 


es ==] Z3 r , ' 
DE DE RE DD TE G9] i 


The axially symmetric potential V (r') which depends on the vector r' and on 
the vibrational coordinates for the nucleus can be expanded as follows: 


VERS VT) V 0(0); (V-4) 
A=0 
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where the Vz (1) are expansion coefficients depending on r and the vibrational co- 
ordinates. It is convenient to introduce a function UQ (7) which represents essentially 
the spherically symmetric part of V (r”), but is slightly modified so as to be inde- 
pendent of the nuclear vibrations. The function V (7) %xr, 44 Yr, Q (&”, p') can then 
be expanded as follows: 


ADDED ENE RED 


1) Fe 1 ( | V-5 
ROEDE DE] VID Dre ROED): | GER 
df Rk 
where 


ve (7) - EN væ (070 ARR 0) 


(Y-6) 
VED 0;0 | GA OV, T)E VU, Y4 79370) | 


according to (V-4) and the fact that the integral (over all directions) of the product 
of three spherical harmonics can be expressed in terms of Clebsch-Gordan coef- 
ficients. If r is larger than the largest axis of the nuclear ellipsoid, the term in (V-6) 
which corresponds to 4 = 2 can be written 


(1 EON 


LEES ULD: OS 0) BEDST De (AVT 


as is easily seen from the formulae (V-10) and (V-13) which will be given later. The 
quantities <44| vr (7) | 4> appearing in (V-5) are identically equal to zero, unless 
GERD) Gr Dis revenfandkbothiltandel farer Q!, and they also have the following 
symmetry properties: 


| VE 202) | 477 = <a | VP (7) | 475 =<ar] V7 2) | 4» 
From (III-3) and (V-5) we get 


VOTE SES HUG) 981.00 


ÅSFel 


Tan lg 7 ' ', T 
Ur Dyk, rO GV, gØ 7) FEB) 


ge VET) (9) IKg, 1, DI 
df v 


140 fol IG ' ' 
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if, in the symmetrization term, use is also made of the fact that <Q | vr? SACD) | qp> is 
independent of the sign of 2 and is zero unless (44+t0)—(944,+1') is even. It is tacitly 
implied that the sum over Ky includes also a sum over the vibrational quantum number 
gf on which REG) depends, although this dependence is not indicated explicitly 
in the notation. 

The potential V (7”) is strongly influenced by the nuclear forces when the alpha 
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particle is inside the nucleus”, but it is due only to the electrostatic forces when the 
alpha particle is outside the surface of the daughter nucleus, and then the potential 
is given by 

PESTO (ERE) 


[r'—r”| dr > (ED 


V(r) -) 
where & is the absolute value of the electron charge, and Q(r'””) is the charge density 
of the daughter nucleus at the point characterized by the position vector 7” which 
(like 7”) is referred to the coordinate system having axes whose directions are fixed 
with respect to the nuclear ellipsoid. From (V-9) we find that, if r is larger than the 
largest axis of the nuclear ellipsoid, the coefficients Vz (1), defined by (V-4), are 
given by 

878 
10) Fe ; SER 
where 
Q3= frå Y, (9) eQ (7) dr". (V-11) 


Ås special cases of this formula, one has 
My SE 3 (V-12) 


since the total charge of the daughter nucleus is (Z—2)8£, and 


ED) ED É 
02 =3|/ ge 0 (V-13) 
where 
8 Qo=Yr? 2 Ps (cos 9") 9 (7”) dr', (V-14) 


Qo being the intrinsic quadrupole moment of the daughter nucleus. 
If the nuclear surface is represented by 


R(W)=R, ( DERE (O)) (V-15) 
4>0 
and the nucleus is assumed to be uniformly charged, one easily finds that 


AED er i GAS 
CE (2?) By B4 GET (V-16) 


provided that the deformation parameters /, (2>0) are very small. In more general 
cases, Qy depends not only on /z, but also on the other parameters fy (X' + 4). 


The wave function Vil, Kr of the daughter nucleus is an eigenfunction of 
the operators Tpart and Trot, the corresponding eigenvalues being Tyart (Kr) and 


& Inside the nucleus, the usefulness of the simple two-body potential V (r”) may be restricted to the 
region near the nuclear surface. 
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Trot (Kp, If) respectively. From (III-1), (III-2), (III-8), (III-9), and (III-10) we there- 
fore get 


(Tyart + Trot) PS | 
mL æn] i Le, 
= € "SIS | Toari (Kr) Uras le ODD ÅR FO DT ERE) 
lås i; (V-17) 
KOTE I; DE, 
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where ? 
i GE 
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For the alpha groups of even-even nuclei (1; = Ky = 0), this formula becomes 


F21](1+1) E. 
Å0,1,2,2 5; > 00, 2 0,7: (V-19) 


From (I11-3), (V-1), (V-2), (V-3), (V-8), and (V-17) we now get the following 
system of coupled differential equations: 


3 re Se 
2 m dr? KEN KE rr vsgE) 
i (V-20) 
FEER LESS; ( j se Å . 
IE » DE <q4 | ver (7) | 947 IK, V, Q (7) +55 År, KØ Ur 1, Q! (7) re O- 
df 1 (øv 
where 
ip hÆLGIEST 
and 
Er, DE Tvart (4) : ((OVEDD) 


In the interior part of the Coulomb barrier, we can in first approximation neglect 
the last term in (V-20), describing the rotational energy of the daughter nucleus. To 
some extent, we can compensate for this approximation by replacing Fr, in the first 
term of (V-20) by Er, Å Ek, where AEr, is a constant which is to be chosen suit- 
ably. Then, (V-20) is replaced by 


EGE É ma 
Forme CAGD)E (Er, FA Ex,) me MØDE D Ver (7) EDR ÆKO, ()=0H 023) 
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if we neglect the vibrations and the possibility for transitions between different bands 
of the daughter nucleus during the alpha particle penetråtion through the potential 
barrier. This system of coupled differential equations is equivalent to the Schrådinger 


equation 
BETTE HV (rn) = Ce AB |w- OB (V-24) 


where 
1 ' 7 , 6 
Dre YE GID) (V-25) 
TE ro £ 
For those alpha groups of even-even nuclei which emerge from the parent 


nucleus in its ground state 0+, and which leave the daughter nucleus in its lowest 
rotational band (I; = Ky= 0), the differential equations (V-20) become 


Re rd er] 49, s, 0 GTS PE nCr) 90, 1,057 0, (V-26) 


NE øl E 
| 2 m dr? BLE) | PEN f 
the functions 94,7 7 (7) being identically equal to zero, unless 2 = 0. 
So far, we have derived the differential equations for the functions GEN eN 


but we need also the differential equations for the functions TRANG; (G)N These latter 


differential equations can be derived either by starting from (III-1) and proceeding 
in a similar way as when deriving the former differential equations, or they can be 
derived from (V-20)F by using (II 9) and" (V-18).Theresultfis 

me n20 É 
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2 m dr 
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where 


Bx, If ET E EN Tvart (K4) og Trot (Gr. I4) . (V-28) 


Obviously, E 
the recoil energy of the daughter nucleus and, hence, no divergency exists between 
Ekpliy and the transition energy Etrang mentioned in the previous chapter. 


KSETG is equal to the sum of the kinetic energy of the alpha particle and 


In the outer region of the r-space where the deviation of V (r') from spherical 
symmetry can be neglected, (V-27) simplifies to 
ENGE 
— UT) SE T) =0. V-29 
2 m dr? UDE, len, LAD ( ) 
For the alpha groups corresponding to transitions from the ground state 0 + of 
an even-even parent nucleus to the lowest rotational band of the daughter nucleus 
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(GE Kr =0OPTET(V-27) becomes 
hak ds EUGEN) 
bane 29 f0,1,1 C) 
FEDE VE TRE DE 
7 


the functions ne (T) being equal to zero, unless /= Is. From (III-11) and (III-12) 
it is immediately seen that the differential equations (V-26) and (V-30) are equivalent. 

The differential equations for the alpha decay of deformed nuclei have also 
been considered by Rasmussen, by RASMUSSEN and SEGALL%, and by STRUTINSKY333, 
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VI. Approximate Solution of the Alpha Penetration Problem. 


In the region just outside the nuclear surface, where appreciable exchange of 
angular momentum takes place between the alpha particle and the daughter nucleus, 
the term TY%ot in the Hamiltonian is relatively unimportant, and we shall neglect it 
in this region”, Neglecting also transitions between different rotational bands of the 
daughter nucleus, we can then in this interior region use the simple three-dimensional 
differential equation?”£ (V-24), the solution of which gives SEE) according to 
(V-25) and ENE (7) according to (III-10). On the other hand, we know that, for 
large values of r the potential V (7”) is nearly spherically symmetric, so that we can 
use the differential equation (V-29) for (Er, (7). In this outer region, the function 
MASSE (T) is therefore approximately equal to a solution of (V-29) representing an 
outgoing wave for large values of r. Such a solution is uniquely determined, except 
for an arbitrary factor which does not depend on r. With a convenient normalization 
of this constant factor, we denote the required solution of (V-29) by G, (Er, t) 


(cf. Appendix A). To obtain an approximate solution of the alpha penetration problem 
we assume that the interior region where (V-24) can be used, and the exterior region 
where (V-29) can be used, overlap. In the common part of these two regions, we 
imagine a spherical surface whose center coincides with that of the nucleus. This 
sphere of radius Ry may be considered as the boundary between the two regions 
where the above mentioned two different approximation procedures apply. The value 
of the radius Ry does not appear in the final formulae, but one should imagine R; 
to be chosen such that the approximations introduced can be justified as well as 


£) For the justification of our approximation procedure it would have been somewhat better to re- 
place Tyot by a suitably chesen constant AEgK, in the interior region instead of neglecting it there. (See 
(V-23) or (V-24) and aiso the footnote on p. 20.) The value of B, defined by (VI-9) and appearing in the final 
formulae, would slightly depend on the value chosen for the constant AEK . This slight dependence is, how- 
ever, unimportant. å 

<£ As has been explained in the previous footnote, the constant AEK appearing in (V-24) is chosen 
equal to zero. 5 
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possible. To this purpose, it seems appropriate to choose Ry such that the above 
mentioned spherical surface lies approximately in the middle of the potential barrier. 
The surface of the daughter nucleus is assumed to be represented byt 


R(&) = Ro da By Yz, 00), (VI-1) 
=0 


where $o=)/47x, $1=0 and |8z|<<|$2]| for 4>2. For small values of $2, the intrinsic 
quadrupole moment Q and the deformation parameter 2 are connected according to 
ihesformula 


3 
Qo = jog (Z—2) Rå B2 90, (VI-2) 


where Z is the atomic number of the parent nucleus and q9 is a dimensionless quantity 
which depends on the charge distribution of the nucleus. The value of 22 can be 
estimated from the formula (VI-2) by assuming the nucleus to have ellipsoidal shape 
and uniform charge density, in which case go is equal to 1. 

The differential equation (V-24) can be written 


Ayp=RKw, (VI-3) 
where" 
V(r')-Ex 
K=K(r')=k É gr: mn (GÆR) 
Kr 
if k here denotes the magnitude of the wave vector corresponding to the energy Ex,. 


According to (V-4), (V-10), (V-12), (VI-4), and (A-4) we get 
1 
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It is convenient to write 
KEDE ER (TEAK (TA (VI-6) 
where 
KØD) Ek Vær (VET 


and AK (Tr, 9”) is the purely anisotropic part of K (7). From (VI-5), (VI-6), and 
(VI-7) we easily find the approximate formula 


alf 
Eg SØEN REDER SE 
ke 2 La 0% Her ae AN MERE , 
AK GES EkE—> D, rl) , z| (z 1) FE GN: (VIEE3) 
10 


+ Our deformation parameter Ø, corresponds to the quantity f cosy in the papers by Boxr and 
MOTTELSON!4> 16 where f is positive and y is equal to 0 for an axially symmetric prolate nucleus and equal 
to ax for an axially symmetric oblate nucleus. 

KRES eet gs VS7)KONEDNSGTinkrert s16 For dEgss (V-13) and (V-16)Finttherpresent” paper: 

ts See footnote Fron pp 16: 
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The differential equation (V1-3), which we use in the region between the nuclear 
surface and the surface of the sphere of radius Ry, is to be solved under the boundary 
conditions that, on the nuclear surface, the wave function wy shall be equal to a given 
function wo (9”, p") depending on Z,A, K;, and Kyp, but not on I; and Ip, and that y 
shall decrease very rapidly as one moves outwards in the potential barrier. An ap- 
proximate solution of this problem can be obtained by arguing as follows. The an- 
isotropic part of the electrostatic field acts essentially only at points comparatively 
near to the nucleus, whereas the centrifugal barrier is of importance over larger distances. 
In first approximation, we may therefore neglect the centrifugal barrier at those small 
distances where the anisotropy is of importance. The fact that this approximation is 
not a very good one is not particularly dangerous, for it turns out that the nuclear 
shape is more important for the aipha decay than is the anisotropy of the electrostatic 
field. From these arguments it follows that, if we want to calculate the alpha wave 
function at a point at a distance r (< R7) from the center of the nucleus, we can in 
first approximation replace the actual nuclear surface by the sphere of radius RY 
and the actual anisotropic potential barrier by the corresponding spherically sym- 
metric potential barrier, if we at the same time replace the alpha wave function 
po (9”, p') on the actual nuclear surface by the alpha wave function given by (B-34) 
on the sphere of radius Ry. In the bracket of (B-34) the first term is due to the non- 
spherical nuclear shape, whereas the second term (which is less important) is due to 
the anisotropy of the sag barrier. For the alpha wave function p at a point” 
(7, &, 9) such that r < Ry we now easily realize that the formula (B-35) applies. This 
formula, which we have here deduced in a more or less intuitive way, is obtained, 
in Appendix B, as an approximate result of a method which, in principle, is more ge- 
neral and less approximate. 

From (V-25), (VI-1), (VI-6), (VI-7), (VI-8), and (B-35) we get 
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Ed 
Here, the unprimed polar angles å and p are assumed to be referred to a coordinate system with 
axes having directions which are fixed with respect to the nucleus. 
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for r < R1. In the exponent of this expression, the term corresponding to 4 = 2 can 
be simplified by means of (V-13) and (VI-2), and for r = R1 the upper integration 
limit kr/x may approximately be replaced by 1. Therefore we get 


. RL (Ekp Ry) Sl (El V 
Ik,,1,2 R1) = 0 Er R,) NL A p) exp 4 f B Po (cos BYRD AYG QW,p)sinddddp, 
g f? 0. 
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where 


and 


Z4>0 
A+2 (VI-10) 
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p (F, 9) = wo (9, p) exp | DSA VS 


the second term in the parenthesis ( ) in the exponent of (VI-10) being less than 20 9/4 
of the first term in this parenthesis if 42>2. The function y, (&, 9) is now assumed 
to be resolved in terms of spherical harmonics as follows 


p; (4, 9) = za, dy, 9 Yy 7 (WP), (VI) 
TØ 
where ay — are coefficients which depend on Er, 2, ANN Den Dagene OF OR OSS sand 


on the parameters appearing in yo (Ø, w). Since the component of the total angular 
momentum along the nuclear axis is conserved during formation of the alpha particle, 
the coefficients a, 9 are different from zero only if 2 = + K; — Ky and | 2] <l. Using 
(III-10), (VI-11), and the formula for IK, 1,2 (R1) on this page, we get 


1x,. 1, Ip (R1) 
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where the quantities 155 v(B) are the elements of a symmetric matrix defined by 


471 
FØCB) = |ø, Q(d) exp iBP2 (cos &)) Oy Q(Y) sin ddd, 
”0 


the normalized Ø-depending factor in the spherical harmonic Y, (9, 9) being denoted 

by 0; 9(9). According to Fig. 6 (in Appendix A), the function G;(E, R1)/G; (E, Ry) 

does not change as rapidly with E in the interior region of the potential barrier, as one 

might expect from the strong E-dependence of G;(E, r). In the functions G, appearing in 
3+ 
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the above formula for fxy,1,174(R1), we may therefore approximately replace" Ex, by Er, 
which is an abbreviation for Ex, 14. Using the formula for fy, 1, 14$(R1) thus obtained and 
noting that forr > Ry the function fxy, 1,14 (7) is approximately equal to G1 (E14, T), except 
for a factor which does not depend on r, we get the formula 


frk4,1,14(T) 
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which is valid for 172 R1. 

From (IV-1), (VI-12), (A-15a), and (A-11b), we find that the alpha transition 
probability per unit time for the transition from the state (K;, I;) of the parent nucleus 
to the state (Kr, If) of the daughter nucleus is 
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According to (IV-2), (VI-12), (A-15a), and (A-11b), the probability per unit 
time and unit solid angle (in a center-of-mass coordinate system with axes having 
directions fixed in space) that an alpha particle leaks out in the direction (2, æ) 
leaving the daughter nucleus in the state (Kyp, I4) is obtained by averaging the fol- 
lowing expression with respect to the initial distribution of M;-values: 
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For K; = 0 the sum over Q in (VI-12), (VI-13), and (VI-14) contains only one 
term (corresponding to 2 =—K7). For K; + 0 the same sum may contain two terms 
(corresponding to Q = K;— Kr and Q=—K;—Ky) if 1> K;+Kp, one term (corre- 
sponding to 2 = K;—K;) if |K;—K/|< 1<K;+ Kyp, whereas it is zero if I<|K:-Ky]. 
Furthermore, / is restricted either to even values only or to odd values only, depending 
on whether the parent nucleus and the band Ky of the daughter nucleus have the 
same or opposite parity. 

A detailed theoretical estimate of the accuracy of the approximations involved 
in the derivation of (VI-12), from which (VI-13) and (VI-14) are obtained, is dif- 


£ The justification for this replacement may be somewhat improved by keeping AEK, in (V-24) and 


choosing it such that EK,” A EKy is some suitable average value of the Elg -values considered. See the foot- 
note + on p. 16. 
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ficult, but a very crude estimate indicates that; for the typical values Z => 90, 
Qo = 10:10 3% cmZ, Ek, æ 5 MeV, and Er, Ek, <= 100 keV, the error involved is less 
than about 50/4, and it seems quite probable that it is appreciably smaller. An 
essential part of this error is due to the approximation that the penetration of the 
alpha particle through the interior region of the potential barrier was assumed to 
take place in the same way as if the nucleus were fixed in space. Another essential 
part of the error is due to the neglect of the anisotropy of the potential barrier in 
the region outside of r = R1. The accuracy of (VI-12) is expected to decrease as 
I4$— K£ increases, since, firstly, the approximation which is involved in the replace- 
ment of (V-20) by (V-24) in the region between the nuclear surface and the sphere 
of radius R1 as well as the replacement of Gi (Ex,, R1) by Gi (E14, R1) becomes less 
accurate as the energy differences between the different alpha groups considered in- 
crease; secondly, in the region r > R1, the neglected coupling between the different 
functions fxæp,1,14(7) should be comparatively more important for the smaller than 
for the larger of these functions. 

In a few special cases, the coupled differential equations which describe the 
alpha decay of even-even nuclei of ellipsoidal shape have been solved numerically 
by RASMUSSEN and SEGALL%, These authors treat the rotational energy term by an 
approximation method and take into account only the partial waves corresponding 
to alpha groups of angular momenta [< 4. Their results may therefore possibly 
contain a minor uncertainty in the intensity ratio for the alpha groups &2 and %&9 and 
a somewhat larger uncertainty in the intensity ratio for the alpha groups &4 and 
%&0. In two cases corresponding to f2-values of 0.30 (their &, = 1.51) and 0.36 (their 
& = 1.41), the alpha wave function was assumed to be constant on the nuclear 
surface. For these two cases the intensity ratios calculated by RASMUSSEN and SEGALL 
(cf. their table I) can be compared with the corresponding intensity ratios calculated 
according to the intensity formula in the present paper, and it is found that the agree- 
ment is good to within 50 ”/,. For the case corresponding to the smaller one of the 
above mentioned f2-values (i. e., to $2 = 0.30), the agreement for the ratio of the in- 
tensities of the alpha groups &2 and &9 is even somewhat better than 10 "/,. 


VII. Analysis of the Experimental Data.” 


The available experimental data on alpha decay of parent nuclei with known 
alpha branching?” and with Z > 82 (except for the data on long range alpha particles) are 
collected in Table 1. Most data are taken from the table of isotopes by HOLLANDER, 
PERLMAN, and SEABORG25, which in the following is referred to as HPS. In those cases, 
where some or all data come from other publications, the last column of Table 1 gives 


+ Note added in proof: Further data are given in the forthcoming review article on alpha decay 
by PERLMAN and RASMUSSEN43å (see their table I). When the analysis in the present paper was made, 
their data had not yet been available. 

+x As regards alpha decay of nuclei with unknown alpha branching, see, e. g., the recent work on 
the light Po—isotopes by ROSENBLUM and TYREN", 
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the appropriate references in terms of figures referring to the list at the end of the 
present paper. HPS is usually not quoted, even if data from this table are used to com- 
plete data taken from other publications. Thus, almost all the half-lives are cited ac- 
cording to HPS. For those parent nuclei which are quoted without any reference, all 
the data have been obtained from HPS. When some remark is required, the last column 
of Table 1 contains a letter which refers to the list of remarks on p. 33. 

In the first column of Table 1, the symbols for the alpha radioactive parent nuclei 
with Z>82 are listed. The third column gives the experimentally measured half-life 
Tæt of each parent nucleus. This half-life may result from different kinds of radio- 
active branching. The second column gives the percentage of radioactivity that is 
alpha decay. For some parent nuclei, column two shows only lower boundaries for 
the percentage of alpha branching. If this lower boundary is rather close to 100, we 
use the mean value of 100 and the boundary in the calculations that are required to 
obtain the values in some of the following columns. If, on the other hand, the lower 
boundary is close to 0, we use the value of the boundary itself. (Example: If the 
second column gives > 90, we have chosen the value 95 (= mean value of 90 and 100). 
lf the second column gives > 10, we have chosen the value 10). When the second 
column is empty, it is understood that the alpha branching is given as 100 9/9. In 
the fourth column, the energies E, (expressed in MeV) of the alpha particles (except 
for the long range alpha particles) are listed. It should be mentioned that the listed 
energies are the alpha particle energies and not the transition energies". This choice 
was made only in order to simplify the numerical work. For most parent nuclei, the 
alpha particle energies are found in HPS or in the references quoted in the last column 
of Table 1; in some cases, however, where these references give, instead, the energies 
of the gamma rays following the alpha decay, the alpha particle energies E, have been 
calculated from the energy EX of the most energetic alpha particles and the excitation 
energies Eexe of the daughter nucleus. This is easily done by means of the formula 


E,= BER — Eexc + E Eexe, 
where M is the mass of the daughter nucleus, and m is the reduced mass of the alpha 
particle and the daughter nucleus. The last term on the right-hand side of this formula 
is sometimes unimportant, and in such cases we have often neglected it entirely. 
For every parent nucleus, the number of alpha particles in a particular group rela- 
tive to the total number of alpha particles is called the relative intensity or the abun- 
dance of the alpha group. The experimentally determined values of these relative inten- 
sities (expressed in per cent) are listed in the fifth column. For those parent nuclei for 
which only one alpha group has been observed this column is empty. For every parent 
nucleus the sum of the percentages in column five should be 100; due to experimental 
uncertainties this is only approximately true for the values listed in Table 1. 


= As will be remembered, the transition energy of an alpha group has been defined as the sum of 
the energy of the alpha particle and the recoil energy of the daughter nucleus. 
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The partial half-life T (expressed in seconds)-of an alpha group is the product 
of three factors. The first factor is the experimentally measured half-life Tio+ in 
column three (expressed in seconds), the second one is 100 divided by the alpha branch- 
ing percentage in column two, and the third one is 100 divided by the relative intensity 
percentage in column five. Column six gives the values of log107', T being expressed in 
seconds”, These values are listed with two decimals, although in general the accuracy 
of the experiments is not as good. Negative logarithms are recognized by an index n, 
so that, e. g., log10T = 8.28, means log10T = 8.28—10. 

If one plots log10T against Ez for the ground state alpha transitions of the 
even-even nuclei, one finds that the points corresponding to a fixed value of the atomic 
number Z of the parent nuclei lie almost on a straight line, except for those parent 
nuclei which are too near to the closed shell with Z = 82 protons and N = 126 neutrons 
(i. e., A= Z+N = 208). This plot is shown in Fig. 1, where the symbol referring to 
each point denotes the corresponding parent nucleus. For those even-even parent 
nuclei (of Po and Em) which have mass numbers A<214, appreciable deviations 
from the straight lines occur, and therefore these nuclei are not included here. Con- 
sidered individually, some of the straight lines may not be very accurately defined, 
but if one considers them collectively and requires that the distances between any 
two neighbouring lines and also the slopes of the lines shall vary regularly with Z, 
all the straight lines are well enough defined for our purpose. In drawing the lines 
one should, of course, take into account the further requirement that those points which 
are most accurately known from experiments must be especially well approximated 
by the straight lines. The straight lines are drawn graphically in such a way thal 
they fulfill as well as possible all these requirements. As is easily seen from the figure, 
the distance between any two successive straight lines decreases slightly and regularly 
as Z increases. The accuracy of the experimental data available does not allow us 
to draw any regular curves that are better than straight lines. The reason why a few 
of the even-even parent nuclei (especially Pu?3?) show more essential deviations from 
the straight lines may be associated with uncertainties in the experimental data for 
these nuclei, and it would therefore be important to re-examine these alpha decays. 

The smooth function of E, and Z which gives the partial half-lives for the ground 
state alpha transitions of even-even nuclei will be denoted by Ty. The analytical 
formula 

logioTo = 4139.8 + 1.83 (Z — 90) + 0.012 (Z — 9072) TIR | CVIr-1) 
— (52,3+0.30 (Z — 90) + 0.001 (Z — 9072), 


where 79 is expressed in seconds and E, in MeV, gives a good approximation of the 
straight lines in Fig. 1. The experimental data on the ground state alpha transitions 
of the even-even nuclei are, in fact, equally well represented by the formula (VII-1) 

+ The values in columns six to nine have been obtained without regard to the signs >, >, and —< 


which appear for some data in columns two to five. 
4 x 
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0.30 0.35 0.40 ; 0.45 0.50 0.55 

E. 2: 

(ed 
Fig. 1. The 10-logarithms of the partial half-lives T (expressed in seconds) for the ground state alpha tran- 
sitions of the even-even nuclei have been plotted against E2 > where lin, is the corresponding alpha particle 
energy expressed in MeV. The chemical symbols referring to the experimental points denote the parent 
nuclei. The straight lines have been drawn graphically such that the points corresponding to a fixed Z-value 
shall be approximated as well as possible by a straight line and that the sequence of straight lines shall be 
spaced as regularly as possible. For the parent nuclei EmZ?"8, Ra2?22, TH226, 77230, 77232, Cf244 and Cf?26 new data 
have appeared after this figure was drawn. The changes in the corresponding experimental points which 
are implied by these new data are small, and therefore we have not re-drawn the figure. For the parent nuclei 
Cf?50, Cf252, and Fm?3, no experimental data were available when this figure was drawn. The points corre- 
sponding to these three parent nuclei have later been included in the figure, where they are recognized 
as small circles. Hence the system of straight lines corresponding to Z = 84, 86,..., 98 as well as the dotted 
line corresponding to Z = 100 (which has been obtained by extrapolating the regular sequence of lines cor- 
responding to Z = 84, 86,..., 98) was drawn before any data on the three parent nuclei Cf25%, Cf?252, and 

Fm?5t were available. 

The signs > and > which appear in column two of Table 1 indicate that, in the above figure, the 
points corresponding to the parent nuclei Pu?32 and Cm?33 may be displaced somewhat downwards (cf. foot- 
note on p. 23). 

After this figure had been prepared for publication, there appeared data on the ground state alpha 
transitions of the parent nuclei Cm?6 (cf. refs. 19a and 43a), Cm248 (cf. ref. 19a), Fm?59% (cf. ref. 43a), and 
Ems (cpFrablelrandref 433)" 
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asS by the straight lines in Fig. 1. The absolute value of the difference between log10 7% 
according to the formula (VII-1) and log10T4 according to the full-drawn straight 


lines in the figure is less than 0.1 if E,” liestbetweent 035Tandt 05 0 ie tik les 
between 4.0 and 8.2 MeV. 
The formula 
1og1070=4140+1.8(Z2'— 90)+0.01(Z — 90) 
— (52.3 +0.3(Z — 90)), 


VE AE 
EG (VIE2) 


which is an approximation of (VII-1), is also useful. For the parent nuclei of the 
isotopes of Th, U, Pu, Cm, Cf, and Fm in Fig. 1, the experimental data on the ground 
state alpha transitions are about equally well represented by (VII-2) as by the straight 
lines in Fig. 1, and for the parent nuclei of the isotopes of Po, Em, and Ra in Fig. 1 
the experimental data are only slightly better represented by the straight lines in 
Fig. 1 than by (VII-2). The absolute value of the difference between log10T'9 according 


to (VII-2) and log1079 according to the straight lines in Fig. 1 is less than 0.2 if VE? 
lies between 0.35 and 0.50, i. e., if E, lies between 4.0 and 8.2 Mev. 

Å priori the smoothed out function Ty is, of course, defined only for even 
Z-values, but it is very useful to define this function also for odd Z-values. Fig, 2 
shows the straight lines obtained for even Z-values in Fig. 1 together with the straight 
lines for odd Z-values found by graphic interpolation. The system of straight lines 
in Fig. 2 defines TQ graphically as a smooth function of E,, and Z both for even and 
odd values of Z. Column seven of Table 1 gives the values of log107'% (with TY expressed 
in seconds) which have been obtained graphically from Fig. 2 by putting Z equal 
to the atomic number of the parent nucleus and E, equal to the alpha particle energy 
of the alpha group considered. These values of log10Tg are listed with two decimals, 
the second decimal being estimated as either 0 or 5. The intrinsic transition prob- 
ability for an alpha group is given by the quotient 79/T which we shall call the 
F-value of the alpha group. From columns six and seven of Table 1 we have calculated 
the values of 10g10(T9/T') to one decimal, and from these values we have obtained 
the F-values which are listed in column eight. The absolute F-values are approximate 
to the same degree as the function 7%, but the relative F-values for the alpha groups 
of any given parent nucleus are much more accurate. For favoured alpha groups 
such relative F-values are very useful. They will be designated as the reduced transition 
probabilities and will be denoted by c and defined as the quotient of the F-value 
for the alpha group considered and the F-value for the special alpha transition of 
the same nucleus for which AI = 0. The F-values appearing in this definition are 
best calculated by means of (VII-1) or (VII-2), for then one avoids errors due to 
the difficulty of reading the graph in Fig. 2 accurately. The reduced transition prob- 
abilities c, obtained by means of (VII-1), for the alpha groups of even-even nuclei 
and the favoured alpha groups of odd-A nuclei are listed in column nine. Ås regards 


26 


PAGE DD 
HAG DG 
5, 


PUNI SL AL SIN I 
MINIS BILLS 
IDI ALI IRI DIL 
ASA GATDIDS SALA 
DITAS 
FALD IASIASDI LAG 
IM ISD SASDASASDSG 
MIG DODGE 
AM AMDSA GADS GA 
BY/Z SD ID SA GIMLE 
UMD GGA 


IDOL ISG) 
MDD OG ÅR 
VA GGIGIA ARR 
TJ DØDØOD 
MOOD OMI] 
WDY/Z 


W/ D) 
e —AØRDS 
"ZODODDD? 
DOG) 
DD) 
en 0.35 040 = at 045 0.50 0.55 
2 


Fig. 2. The straight lines corresponding to even values of the atomic number Z of the parent nucleus are 
the same as those in Fig. 1; the straight lines corresponding to odd values of Z have been obtained by 
graphic interpolation between the lines corresponding to even Z-values. 


the numerical values of c compared with those of F it is clear that, on the average, 
F x c for the alpha groups of even-even nuclei and, as shown in Table 1, F = 0.5c 
for the favoured alpha groups of odd-A nuclei. 


Note added in proof: In the forthcoming review article on alpha decay by PERLMAN and Ras- 
MUSSEN43a, an analysis of the empirical alpha decay data is also presented. It is essentially similar to 
ours, although it is partly based on quite recent data which had not been available to us. Despite some 
minor differences between the two analyses, the results are in general in satisfactory agreement. 
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TABLE 1. Analysis of the experimental data on alpha decay for parent nuclei 


with Z>82 (except for the data on long range alpha particles). 


HE ERE SERRA SS FRS BEES be 
ssBi <198 7 me 2.01 0.95 0.08 
Bilss 5:10—? 7 m| 5.83 5.92 2.60 0.0005 
Biles 10— Ron m| 5.47 ag 4.30 0.001 
Bi20 3.10—3| 62 wa) ls 8.09 5.95 0.008 
Bizo3 | 1075 12 ne EA .85 11.63 7.65 0.0001 
Bi2os SOE 8855 24.98 21.55 0.0004 
Bi2u 99.68 2.16 m| 6.618 | 84 2.19 9.20, | 0.001 
627250 HET 2.91 0.65 0.005 
Bi21? 389 60.5 mu 6.086 FR 4.60 1.50 0.0008 
6.047 | 69.9 4.19 1.60 0.003 
5.765 ill 5.80 2.85 0.001 
5.622 0.15 6.86 3.55 0.0005 
5.603 er 5.99 3.65 0.005 
5.481 0.016 7.83 4.30 0.0003 
Bi213 2 46.5 m| 5.86 5.14 2.45 0.002 
Bird 0.04 19.7 m| 5.505 | 45 6.82 4.15 0.002 
5.444 | 55 6.73 4.40 0.005 
ssPo?u |x 1 888 mi ESer 6.13 5.30 0.2 
Po2os |x 10 9 ik. BNG 6.89 6.10 0.2 
Pork |O0= 1 05 7 je BO 8.31 6.80 0.03 
Po208 | 2.93 y | 5.109 7.97 6.70 0.05 
Po? |> 90 1100 fk | FEST 9.52 8.05 0.03 
Po | | 139 d "5.301 7.08 5.60 | 0.03 
Po?n 0.52 st 7.434 99 9.72, 6.85, | 0.001 
6.88 0.50 2.02 8.80, | 0.0006 
6.56 0.53 1.99 0.00 0.01 
Po»? BI0STOREESE R8. 776 3.487 0, | Om 
Po213 eos 85336 4.62, 43005 
Po2u4 5 057 680 GELSE Ol, | dl 
Pos |—> 99 1.83-1072 s | 7.365 HØ. HUD | (6 
Pos |—> 99 0.158 SA E GE 9.20, VD. | dl 
Pos |—> 99 3.05 m|… 5.998 2.26 2.25 1 
FREE SET 0 RE 0 HU ES 75 4.86 3.80 0.08 
At20s ere 407 h| 5.65 6.09 4.25 0.02 
eee Er 5.5 25:65 5.60 4.25 0.05 
At210 ig 8.3 ie I EO | SP TATA 4.90 0.002 
5.437 | 31 7.75 5.30 0.004 
5.355 | 37 7.68 5.80 0.01 
At2u 40.9 ÅL h | 5.862 4.81 3.30 0.03 
Atat 202510 ERE SE HR EL TE 4.30, 385, od 
At215 10 SN ES-00 6.00, 5.50, | 0.3 
At21s 5510 HEER E7.79 6.48, Gore Ron 
At217? 0.019 $ | On SE 305 | 25 
At | 97 0.9 m| 6.27 1.75 1.50 0.6 
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Parent 
nucleus 


208 
ssem 
Em?09 
Em?10 
Em?U 


Em2?1? 
Em 
Em?2?16 
Em2?17 
Em 


Em”? 


EmZ?20 
Ems 


Em??2 


7 By IE 


Fr?18 
Fr?19 
F"r?20 


Fr??1 


ssRa213 
Ra219 
Ra220 
ktasel 
Ra222 


Ra?22 


Ra? 


Ra?26 


ol, x 


x 20 
» 20 
> 95 

AD 


44 


w 107 
SIDES 
30 


m 


> 


an 


d 


d 


E& 


Rel. &-int. 
(%/,) 


33 
67 


99.8 


37 
39 


94.3 


log, T 


3.84 
3.95 
3.96 
5.84 
5.54 
3.14 
4.007 
6.00, 
7.007 
85580 
0.98 
0.75 
1.42 
1.51 
1.99 
1.74 
3.87 
5.52 


3.85 
3183 
4.04 
TET OR 
830 
1.44 
2.59 
3107 


DT 
7.00, 
8.487 
1.48 
1.60 
2.88 
7.03 
6.26 
6.60 
7.03 
7.68 
7.51 
5.52 
6.80 
10.73 
11.95 


c 


0.34 


1.35 


References 
and remarks 


11.42 
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TABLE 1 (continuéd). 
ERE Trot SAD ER KEE SO ERR RE 
s9Åc?2? 5.5 S 6.96 0.74 0.35 0.4 
Ac223 99 29 m| 6.64 212 1.60 0.3 
Ac?4 |. 10 2.9 h 6.17 5.02 3.55 0.03 
Ac22?5 10.0 d | 5.80 5.94 5.35 0.3 
Ac227 f£2 22 y 4.942 10.76 10.05 0.2 
s0Th223 01 S 7.55 9.007, SØS] 105 
Th224 gl s 113 0.00 0.10 i 
Th??53 |. 90 8.0 m| 6.57 DENS 2,30 0.4 
Th225 30.9 m (6:33 79 3.37 3.30 0.8 1 le ØE Te 
6.22 19 3.99 3.80 0.6 0.74 
6.095+ 1eg 5.03 4.40 0.3 0.25 
6.03 0.6 5.49 4.65 012 0.17 
Th?27 18.6 da KG: 030 19 6.93 4.70 0.006 c 
6.001 5 BEST 4.80 0.002 
5.972 21 6.88 4.95 0.01 
5.952 13 7.09 5.05 0.01 
58022 SED 7.90 5.20 0.002 
| 5.860 4 7.60 5.50 0.008 
| 5.796 2 7.90 5.85 0.01 
| 5749817 6.97 6.00 0.1 
| STS DR og 6.15 0.008 
| 5.704 15 7.03 6.25 0.2 
DØBTE n | 259) 7.90 6.55 0.05 
Th?28 1.90 y 5.421 Th 7.93 7.80 0.8 1 AES SER 
5.338 28 8233 8.25 0.8 1.15 
5.208+ 0.4 10.18 9.00 0.06 0.089 
5.173 0.2 10.48 9.20 0.05 0.074 
5.137+ 0.03 11850 9.40 0.01 0.018 
rn ze | 7340 y 5:02 —10 12.36 KOS KO) 0.005 
| AOA ER 399 12.06 10.60 0.03 
ASB 70 11.52 005 
Th230 8104 y 4.682 75.5 12.52 12.30 0.6 1 54, b 
1.614 24,5 13.01 12.80 0.6 0.98 
AGO 015 15.23 13.80 0.04 0.062 
4.427% | 0.03 15.92 14.15 0.02 0.028 
Th?32 1.39-1019 y)| 3.98 17.64 17.80 1.6 
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TABLE 1 (continued). 
SEE EEEEEEEELEÆEELL]EYVST STE er FOO AOA TS 


RE: References 
Ses 91, & Ttot E« ETS log, T 10810 T, EK Cc and remarks 
Eee TS 
51.225 1.8 m 6.81 2.03 1:70 0.5 
Pa??? |x 85 38.3 m 6.46 3.43 3.20 0.6 
kbazze »2 22 h 6.09 75 6.72 4.80 0.01 
5.85 25 7.20 6.00 0.06 
Pa?31 3.43:10% —y 5.0490 12 12.95 10.55 0.004 31, a 
5.0205 23 12.67 10.65 0.01 
5.006 26 12.62 10.70 0.01 
4.974 1.5 13.86 10.95 0.001 
4.942 24 12.65 iglgks 0.03 
4.925 3 13.56 11.25 0.005 
4.848 ir5 13.86 11875 0.008 
4.727 10 13.03 12.60 0.4 1 
4.704 0.8 14.13 12.70 0.04 
4.671 1.3 13.92 12.95 0.1 0.32 
4.630 0.3 14.56 13.25 0.05 
ØDE 13 m 6.8 1.89 2.20 2.0 
WER » 80 9.3 m 6.67 2.84 2.70 0.8 
LUE » 20 58 m 6.42 4.24 3319 0.3 
U230 20.8 d 5.884 6751 6.43 6.30 0.8 1 11,42,'be 
5.813 32.1 6375 6.65 0.8 1 É 
5.658% 0.4 8.65 7:90 0.08 0.091 
5.658 0.4 8.65 7.50 0.08 0.091 
Wy23r DES SDER 4.3 d 5.45 9.83 8.60 0.06 
uz 74 y 5.318 68 9.54 9.35 0.6 1 NORET OED 
5.261 32 9.86 9.75 0.8 150 
5.134 0.32 11.86 10.50 0.04 0.059 
4290850 00] 13:33 1835 0.01 0.013 
uæ 1620 REE 4.8157 | 83.5 12.79 152555 0.6 1 31 
| HR 77328 EM Æ9 13.53 12.80 0.2 0.34 
| 1042/1746 14.50 13.20 0.05 0.093 
| 4.4890 | 0.03 14.23 14.85 0.04 
use 2.5:105 y 4.763 LD 13.02 12.95 0.8 1 SÅ EgED: 
4.710 25 13.50 13.25 0.6 0.74 
4.593 0.3 15.42 14.10 0.05 0.062 
25 Ta SK OD y 4.58 10 15835 14.20 0.0008 
4.47 « å 17.87 15.00 0.001 
4.40 83 16.43 15:50 0.1 
4.20 4 | 17.75 187815 0.3 
(225 2.4107 y 4.50 14.88 14.85 1 
(BES AGE LORS EAST 8 78 17.26 17.30 1 het | 54 
| 4.13 22 ES 720 0.8 0.80 
4.02 « 0.1 20.15 18.70 0.03 0.033 
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TABLE 1 (continued). 


ere er Bee 
| 
s3Np?81 550 m| 6.28 3.48 4.85. | 24 
Npåss (1072 | 35 m| 5.53 8.32 8.65 2.0 
Np2ss | 3.5-10—3) 410 da 06 12.00 11.50 0.3 34 
Np237 | 2.20-106 y| 4.872 3 15.35 12.70 0.003 38 
4.816 3.5 15.30 memo 0.006 
ARS 7 ESS 14.12 13.25 om 
4.767 | 29 14.38 13.40 0.1 
4.713 kr, 15.61 13.80 0.02 
4.674 313 15.32 14.05 0.05 
4.644 6.0 15.06 14.30 0.2 
4.589 0.5 16.14 14.70 0.04 
4.52 0.02 17.54 15.20 0.005 
ssPu?3? | > 2 36 m| 6.58 5.03 4.00 0.1 
Pus: | — 4 æg h| 6.19 5.91 5.75 0.6 
Puzss | x 0.002 | 26 m| 5.85 7.89 7.45 0.4 
Pu?3s | DES ill | TE 7.93 7.95 1 
Pu228 | 90 ONE 54050 72 9.60 9.40 0.6 1 TS 
| 5.452 | 28 10.01 9.70 0.5 0.69 
15352 CO ORE 255 08 REN 6258 010 06 RR 000 8g7 
Pu?39 24360 75150 69 127050 11550 0:3 i Sl an 
513700 HE20 12.59 11.55 0.1 0.33 
| 5100900 HEAT 12.84 11.80 0.1 0.35 
Pu240 | 6300 mi 62 76 11.42 11.35 0.8 1 24, 31, 35, 50 
SR | 527 11.92 170 0.6 0.62 
5.014 0.1 14.30 12.35 0.01 0.012 
Pu?u | — 5-10723| 14 ve HE s0s 75 18707 13.15 leR 50, g 
ARS ASER FE25 13.55 13.50 1 
Pu2?d? 5-10: 1 E4:308 0 30 13.30 13.15 0.8 1 50, h 
4.854 |. 20 13.90 13.45 0.4 0.50 
ssAm?37 (005 ES 73 h | 6.01 7.97 7.10 0.1 
Am?39 00037 72 05.75 9.16 8.45 0.2 
Am24l (1470 y| 5.5408 | 0.39 12.58 9.60 0.001 6, 31, i 
| 5:5082 | 0.24 12.79 9.85 0.001 
5.4820 | 85 10.24 9.95 0.5 1 
5143000 1284 11.06 10.25 0.2 0.31 
5.3860 | 1.66 11.95 10.40 0.03 0.083 
5.3210 | 0.015 14.00 10.95 0.001 0.0018 
5.241 0.002 14.87 11.40 0.0003 | 0.00076 
Am?2a3 7.6-1082 y)| 5.340 0.17 14.15 10.85 0.0005 6, 51, j 
5.309 0.16 14.18 11.00 0.0006 
526700 87.1 11.44 11.30 0.8 1 
502400 ETS 12.32 11.55 0.2 0.25 
5.169 ift 18854 11.90 0.04 0.052 
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TABLE 1 (continued). 


ARR | | References 
SR 0, Ttot | Ex Sega log» T | 10870 T, TEN | Cc | and remarks 
EEN ERE ES SE ER 
sgCm?238 | > 10 2.5 h| 6.50 4.95 5.20 1.6 
Cm240 26.8 al |» TØS 6.36 6.35 1 50 
Cm241 0.2 35 Gl |. BERO 9.18 8.20 0.1 50 
Cm?24? 162 d | 6.110 TET 128 MED 0.6 1 3, 50, b 
6.066 26.3 7.73 7.30 0.4 0.60 
5.965 0.035 10.60 7.85 0.002 0.0028 
Cm233 35 y| 5.985 6 10.26 70 0.003 3, 42, 50,08 
5.777 81 9.13 8.85 0.5 1 
5.732 13 9.93 9.10 0.2 0.29 
51675 | Kæ2 10.74 9.40 0.05 
Cm: 19 y 5.798 75 8.90 8.75 0.8 1 3, 50 
5.755 25 9.38 8.90 0.3 0.56 
Cm245 26:103 Un H5:32 11.91 11885 0.3 50, 1 
o7Bk2ss | 0.1 4.6 be 62 30 7.74 4.65 0.0008 
| 6.55 53 7.49 5.45 0.01 
| 6.20 17 7.99 7.05 0.1 
Bk245 0.105 4.95 da WEGS7 33 9.09 6.30 0.002 39 
| REG 7 RET 9.00 7.30 0.02 
(ff 5-89 FN 226 9.19 | 78.60 0.3 
Bk29 | 0.001 | 290 dd 5-40] | 12.40 11.50 0.1 37 
28Cf244 | — 25 mL? Sus ES 50 1 42 
Cf245 mer mi Byer 3.42 3.45 1 42 
Cf246 3587 fa 65758 78 5.22 5.05 0.6 1 36 
| an IE sø BE 5.20 0.3 0.44 
636087 0308 SER 5.65 0.003 0.0046 
Cf238 205 d| 6.26 7.25 7.30 1:3 50 
Cf239 500 ul ET 3 HE E72 7.60 0.00008 29 
5:08 MEG (11.42500 0:00 0.004 | 
5.80 | 90 (HET O24 RE EOS 70 0.3 | 
Cf250 SE ul 6024 | "83 18853 8555 1 | 10 
| 5.980 17 | HE50'27 8.70 0.3 0.35 
Cf252 KELD UN 6-112 84.5 (97:87 8.05 1.6 1 10 
6.069 15.5 8.61 8.25 0.4 0.30 
99253 | 19.3 d | 6.633 90.2 6.27 6.00 0.5 1 42 
| 6.592 TET (8:34 6.10 0.06 0.13 
6.545 1.8 TRO 6.30 0.02 0.050 
6.493 0.3 REE 2:74 6.60 0.008 
6.25 0.05-0.1| 9.35 7.80 0.03 
100Fm?82 DL HS 70450 FEET reg | 0.5 | 26 
Em?2s4 3.3 På 0717000 85 (14 105 0.8 1 (er hr 589S 
HELDS 15 | 4.90 4.15 0.2 0.26 | 
| 7.034 |>0.2 6.77 4.60 0.006 0.0085 | 
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REMARKS TO TABLE 1. 


BRIGGS!” gives quite different data for RaZ?23, 

The odd parity alpha groups of even-even nuclei are recognized by asterisks + placed after their E&- 
values. The assignment 1— for the lowest odd parity state of the daughter nucleus is certain for the 
parent nuclei Th?26, Th228, TH230, 17230, and almost sure for the parent nucleus U?22, In the alpha decay 
of Cm??, one has also found a 0+ — 1— alpha transition (not listed here), but, in this case, the 1— 
interpretation is somewhat doubtful??, See also Note added in proof, p. 46. 

The data on Th?” reported by FRILLEY et al.?? are the same as those on Th??? in HPS which are listed 
first, but BRIGGS”! data are different. 

The two alpha groups of energies 4.925 MeV and 4.630 MeV were not reported by GoLDin et al.3; they 
are obtained by HUMMEL (private communication from Professor PERLMAN). 

The two alpha groups of 5.658 MeV energy have actually not been resolved, but beta decay data indicate 
that about half of the 5.658 MeV alpha particles leave the daughter nucleus in the 4+ state and the 
other half in the 1— state (private communication from Professor PERLMAN). 

The most energetic alpha group which has been observed for the parent nucleus Pu?29 (i, e., the group 
with the alpha particle energy 5.150 MeV) does not appear to correspond to a transition to the ground 
state of the daughter nucleus?. 

The decay energy of Pu?" (calculated from closed cycles) is 5.13 MeV, according to GLass et al.2?, The 
most energetic alpha group which has been observed for the parent nucleus Pu?% (i. e., the group with 
the alpha particle energy 4.893 MeV) has a transition energy which is less than the decay energy by 
about 0.150 MeV, and therefore it seems probable that this alpha group does not leave the daughter 
nucleus in its ground state. 

For the half-life of Pu?”2?, we have used the value 5-10 years, which is given in HPS. The value 9: 105 
years given by SEABORG and Katz? js too high. New measurements have yilded a value of about 
45102 years" 29 22; 

It should also be remarked that the measurement (by AsaRro) from which the values 80 ?/, and 
20'/, for the relative intensities of the two alpha groups of the parent nucleus Pu?? are obtained is 
very rough (private communication from Professor PERLMAN). 

There may be some uncertainty about the 5.241 MeV alpha group. 

On the basis of experimental indication for a ground state spin 1/2 for Np?3%, objection has been raised 
by STEPHENS et al. against the assignment of the spin 5/2 to the level in which the alpha group with the 
alpha particle energy 5.267 MeV leaves the daughter nucleus Np?23?, The assignment of spin 1/2 for the 
ground state of Np?3% implies, however, many difficulties for the interpretation of the level structure 
of Np?3? "as well as the beta decay of this nucleus 25". 

The 5.985 MeV alpha group from the parent nucleus Cm?” leaves the daughter nucleus in an excited 
state which lies 76 keV above the ground state %2", 

The state in which the 5.675 MeV alpha group leaves the daughter nucleus decays to the lowest 

state of the favoured band (in which the 5.777 MeV alpha group leaves the daughter nucleus) by an 
E1 gamma transition, and therefore the first of these states probably consists of two different states 
lying close together (private communication from Professor PERLMAN). This fact modifies an inter- 
pretation which has recently been given by A. BouRr et al.!5, 
For the parent nucleus Cm?%5, the data given by SEABORG and Kartz%? are more recent and probably 
more accurate than those given in HPS. It should also be remarked that alpha-gamma coincidence 
experiments indicate that the 5.34 MeV alpha group does not leave the daughter nucleus in its ground 
state (private communication from Professor PERLMAN). 


Mat. Fys.Skr. Dan, Vid, Selsk. 1, no.3. 5 


34 


ABE DE 


Ness 


For the alpha groups of even-even nuclei, which are characterized by their an- 


gular momenta l, we here list the empirical hindrance factors which we have obtained 
as private communication from Professor PERLMAN. The reduced transition probabi- 
lities cz have been calculated as the reciprocals of these hindrance factors (except for 
the alpha group &3 of the parent nucleus Th?28 for which cz has been obtained from 
Table 1). From the c;-values we have then calculated the b;-values according to (VIII-7) 


and (VIII-12). 


| 
| 


Parent Hindrance Parent Hindrance 
nucleus factor cl bi nucleus factor dl bi 
ange 2 1.20 0.83 IETES vs P) 15 0.67 1.06 
Ra222 2 1220 0.83 LØ suse » 125 0.67 1.05 
Ras PSR 0.90 ll 1.38 [yes 2 1.6 0.62 102 
| 4 115 0.0087 0.22 
she. PA 1.50 0.67 1.06 6 385 1 0.0026 0.30 
1 4.0 Of25 0.55 8 15000 0.000067' 0.18 
Å 5 O13 0.84 
Puz 2 127 0.59 0.99 
une 2 0.90 tsk 153% 4 68 0.015 0.28 
1 10 0.10 0.35 | 
4 | 12 0.083 0.69 CNE 2 1.65 0.61 1.00 
3 0.018 0:23 4 380 0.0026 O2 
6 260 0.0038 0.36 
oe 2 1.0 1.00 1230 8 4850 0.00021 0.28 
4 9.8 0.10 (MED AFÅ i: 1000 0.0010 0.035 
1 40 0.025 01% 
6 |> 8500 0.00012 < 0.070 Cm D 1.8 0.56 0.96 
4 1000 0.0010 0.073 
inn 2, 0.80 125 1.47 6 340 0.0029 0.31 
SØDEED 2 il 0.91 1523 orll22s 2 2.7 OS 0.77 
1 11 0.091 0.33 4 115 0.0087 0.21 
4 9.6 0.10 0.76 6 320 0.0031 | 0.31 
bær » et 0.91 1294 Cf250 D) 380 0.33 03 
4 16 0.062 0.60 
ile] 70 0.014 0.13 Ci 2 32 URE OSSE 0.71 
& | 270 0.0037 0.37 4 — 94 — 0.011 — 0.23 
U?34 & 162, 0.83 1618 EN RU UGE 52 4.0 HO 25 | 0.63 
4) 14,5 0.069 0.62 4 54 | 0.018 0.30 
| 6 > 800 < 0.0012 Sa UEr K0) 
mee 2 | lg! 0.91 1.24 


=£ Note added in proof: Somewhat more recent data are given in the review article on alpha decay 
by PERLMAN and RASMUSSEN?" (see their fig. 15) which will appear in Handbuch der Physik. To our 
purpose, however, the data in Table 2 can be used equally well. i 
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In Table 2 we list empirical values of the hindrance factors for the alpha groups 
of even-even nuclei, kindly communicated to us by Professor PERLMAN. This table 
contains some hindrance factors for alpha groups populating higher states of the 
daughter nucleus which cannot be obtained from Table 1. The reduced transition 
probabilities c; in Table 2 have been calculated as the reciprocals of the hindrance 
factors listed there (except for the alpha group &3 of Th?28 for which c3 has been ob- 
tained from Table 1). With a few exceptions, the c;-values in Table 2 are in satis- 
factory agreement with those in Table 1. The exceptions may be associated with 
differences in the definitions used for the hindrance factors and their reciprocals, 
which are the reduced transition probabilities c;. For the discussion of the empirical 
data in relation to the present theory, which will be given in the next chapter, it is 
better to use not the reduced transition probabilities c; themselves, but certain quan- 
tities b; which are defined by (VIII-12). The b;-values which are listed in Table 2 
have been calculated from the c;-values there by assuming the nuclear radii to be 
svend kby (MILL 7): 


VIII. Comparison between Theory and Empirical Data. 


In this chapter, we first give a brief review of the theoretical treatment devel- 
oped in the Chapters III—VI, and subsequently discuss the ability of this theory to 
account for the available empirical data. 

To treat the alpha fine structure problem we must describe the penetration of 
the alpha particle through the potential barrier around the deformed daughter nucleus 
by means of a wave function which takes into account the possibility for the alpha 
particle to leak out in different channels corresponding to different final states of the 
daughter nucleus. Using the center-of-mass coordinate system for the alpha particle 
and the daughter nucleus we have constructed this wave function in such a way that 
its total angular momentum is I; and its angular momentum component along 
a fixed direction in space is M;, where I; and M; denote the corresponding quantum 
numbers for the parent nucleus. This wave function contains for each channel a 
function which depends on the distance r of the alpha particle from the center of the 
daughter nucleus and which will be referred to as the radial channel function. These 
channel functions satisfy a system of coupled differential equations, which has been 
derived on the assumption that the daughter nucleus possesses axial symmetry 
and can be described in terms of a rotational motion superimposed on the intrinsic 
motion of the nucleons in the deformed nuclear potential. If the small coupling causing 
transitions between different rotational bands of the daughter nucleus during the alpha 
particle penetration through the potential barrier is neglected, the differential equations 
couple only those radial channel functions which correspond to alpha groups which 
emerge from a given state (K;, 1I;) of the parent nucleus and leave the daughter nucleus 
in the same rotational band characterized by the quantum number Kyp. In terms of 
the alpha wave function on the nuclear surface, these differential equations determine 
the transition probabilities for the different alpha groups. 
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The passage of the alpha particle through the region near the nuclear surface, 
where appreciable exchange of angular momentum occurs between the alpha particle 
and the daughter nucleus, takes place in a time which is short compared to the nuclear 
rotational period. Therefore it is in first approximation allowed to consider the nucleus 
as fixed in space, while the alpha particle traverses this interior region. An approximate 
solution of the problem concerning the penetration of an alpha particle through such 
a potential barrier fixed in space can be obtained by means of an approximation 
procedure which has been described in a qualitative way in Chapter VI, and which 
is developed more rigorously in Appendix B. Using this approximation method, one 
finds that in the interior part of the potential barrier of a deformed nucleus the alpha 
wave function can in first approximation be calculated from the alpha wave function 
on the nuclear surface by means of a formula which represents a generalization of 
the well-known formula for the spherically symmetric potential barrier. 

In the outer part of the potential barrier and outside of this barrier the Coulomb 
field around the daughter nucleus is almost spherically symmetric, and there is then 
practically no coupling between the different channels into which the alpha particle 
is able to leak out. Therefore, the passage of the alpha particle through this outer 
region takes place approximately in the same way as the passage of an alpha particle 
of given energy through a spherically symmetric potential barrier, the alpha transition 
energy being that corresponding to the channel considered. Except for constant factors, 
one can therefore in the outer region easily find the radial channel functions which 
correspond to outgoing waves. 

Assuming that there exists an intermediate region around a spherical surface of 
radius Ry, where the above approximations for the interior and exterior regions both 
apply, we can determine the unknown constant factors which appear in our approx- 
imate solution for the outer region by matching this solution to that for the interior 
region at the spherical surface r = R. In this way, the radial channel functions are 
determined in the whole space outside the nuclear surface. They depend, of course, 
on the boundary conditions for the alpha wave function on the nuclear surface. 

By applying the approximation procedure outlined above to the alpha transition 
from the state (I;, K;) of a given parent nucleus to the state (Ip, Kr) of its daughter 
nucleus, we have in Chapter VI found that, for r > Ry, the radial channel function 
fx4,1,14(T) corresponding to the angular momentum I is given by (VI-12), the 
transition probability Px,,174 is given by (VI-13), and the angular distribution of 
the alpha group considered is given by (VI-14). Thus it is seen that the alpha decay 
depends on a quantity B, defined by (VI-9), which is a measure of the quadrupole 
deformation, and on the expansion coefficients a; of the function p1, defined by (VI-10), 
which depends on the alpha wave function on the nuclear surface, 9, as well as on 
the higher multipoles of the nuclear surface. The further dependence on the transition 
energy, the atomic number, the nuclear radius, and the quantum numbers Ky, I;, 
and I£ should also be mentioned. 

Instead of considering the transition probabilities it is convenient to introduce 
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the so-called reduced transition probabilities in which the strong energy dependence 
has been separated off. In the definition used in Chapter VII, this energy dependence 
has been taken directly from the empirical data on the ground state alpha transitions 
of even-even nuclei. This definition of the reduced transition probabilities is con- 
venient when one deals with the empirical data, but in more theoretical connections it 
is often advantageous to choose a slightly different definition according to which the 
energy depending factor to be separated off is taken from the well-known theory for the 
penetration of an alpha particle of angular momentum zero through a spherically sym- 
metric potential barrier. For practical purposes, it is in general unnecessary to distin- 
guish between these two definitions, for they give approximately the same numerical 
results, unless the energy differences of the alpha groups considered are too large. 

In the following sections, we chiefly discuss the alpha decays of even-even nuclei 
and the favoured alpha decays of odd-Å nuclei. For these alpha decays the alpha 
particle is formed with Q = 0. In order to simplify the notations we shall, therefore, 
in general write k, ,(B) instead of k? ,,(B) and a, instead of da, 0: 


a. Ground state alpha transition in even-even nuclei. 


According to Chapter VI, the alpha decay of an actual nucleus takes place in 
the same way as that of a fictive nucleus with the purely ellipsoidal surface 


BOV ERE (IR Y20(2)) (VIII-1) 


on which the alpha wave function is equal to y1(9") and outside of which the electro- 
static potential is assumed to be equal to the sum of the spherically symmetric poten- 
tial and the pure quadrupole potential of the actual nucleus. The angular dependence 
of p1 (9) is due both to the higher multipoles in the nuclear deformation and to the 
variation of the alpha wave function on the nuclear surface. For the ground state 
alpha transition probability of an even-even nucleus we obtain 
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as a special case of (VI-13). The quantities a; appearing here are the coefficients in 
the expansion of p1 (9”) in terms of the spherical harmonics Y7,9(9'). From (VIII-2) 
we find that, if we assume y1 (9”) to be independent of % and to remain constant 
while 2 varies, the ellipsoidal deformation of the nucleus causes the ground state 
alpha transition probability to be enhanced by the factor” 
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£ The enhancement factor (VIII-3) has been obtained on the assumption that R, remains constant 
as P, varies. There is a slight difference between this enhancement factor and that corresponding to constant 
nuclear volume for varying $,, but for our purpose the difference is unimportant. 
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which increases as |B| increases, 1. €., as | 62 | increases. HInL and WHEELERZ3, who 
were the first to point at the importance of the nuclear deformation for alpha decay, 
estimated the enhancement very roughly and obtained an expression which in our 


notations reads 
Max era Pieor De: (VIII-4) 
ø 


however, this expression overestimates the enhancement considerably". By means of 
a straightforward improvement of the considerations used by Hizx and WHEELER 
one may replace their expression for the enhancement by 


O 1 
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É 47 
but also this expression overestimates the effect, as is easily seen by comparing (VIII-5) 
with (VIII-3) and using ScHwaRrz's inequality. 
For difterent values of B, we give in Table 3 the deformation parameter 2, the. 
intrinsic quadrupole moment Qg, and the value of the enhancement factor according 


TABLE 3. Calculated values of the enhancement factor for different 
nuclear deformations and in different degrees of approximation. 


Enhancement factor 
ME Ve 8. B according to 


(VIII-3) (IFs Hr 

88 — 0.35 mz"3 3.49 5.93 20.1 

1985 HE 0:94 73) 1.88 SY 7.4 

=63 = ØD æg 1.20 37 DT 
0 0 0 1 1 | 1 

6.3 0.12 1 1.24 1.55 | 7.4 

12.5 0.24 2 2.41 5.11 54.6 

18.8 0.35 3 6.93 23.98 403.4 


to the correct formula (VIII-3) as well as according to (VIII-4) and (VIII-5). The 
values of $» have been obtained from the approximate formula 


B = 8.5 f> (VIII-6) 


which is found to be fulfilled to within 10 /,. This formula follows from (VI-9) and 
the empirical alpha decay data for even-even nuclei, if the nucleus is assumed to be 


£ On the basis of numerical calculations, the increase in the alpha decay rate due to the nuclear 
deformation has been discussed by RASMUSSEN and SEGALL?? who obtained much smaller enhancement 
factors than Hix1 and WHEELER%, Their results seem to be consistent with ours, but a detailed comparison 
is difficult because of different boundary conditions for the alpha wave function. 
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uniformly charged (i. e., go = 1) and the radius Ry of the daughter is assumed to be 
given by 1 
Ro l44E 0 (44) erne (VIII-7) 


where Å is the mass number of the parent nucleus. The values of the intrinsic qua- 
drupole moment Q% given in Table 3 have been calculated from (VI-2) and (VIII-7) 
by assuming that g0 = 1, Z—2 = 90, and A—4 = 230. 

There is a certain ambiguity in the choice of the radius Ro, and there may be 
some reasons for choosing Ro somewhat smaller than according to (VIII-7), e. g., 


till 
hy SKO PS IND orne (VIII-8) 


The choice of such a smaller radius would, however, not change the qualitative con- 
clusions which we shall arrive at in this chapter. The larger radius (VIII-7) appears 
to represent the extension of the range of interaction between the alpha particle and 
the daughter nucleus?”, For the estimate of 2 from the electric quadrupole moment 
Qo a somewhat smaller radius should be employed. With the present accuracy of 
the Qo-values this point is, however, of minor significance. 

Empirically it is known that the probabilities for the ground state alpha tran- 
sitions of even-even nuclei can be represented with good approximation (within a 
factor of about 2) by a smoothed out function Po (Z, EQ) which depends only on the 
atomic number Z of the parent nucleus and on the alpha decay energy EQ. This 
function can easily be obtained from the formula 
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(VIII-9) 
where Tp is the smoothed out half-life function (depending on the atomic number Z 
of the parent nucleus and on the alpha particle energy E4) which has been determined 
in a semi-empirical manner in Chapter VII. The variation of the enhancement factor 
with the nuclear deformation should obviously affect the function Po (Z, EQ). From 
Coulomb excitation data it has been found that Qo > 10 for the heavy elements around 
U. (See the review article by ALDER et al.7.) For these QQ-values the correct enhance- 
ment factor (calculated from (VIII-3)) is about 2 according to Table 3. A factor in 
Po (Z, EQ) of at most this order of magnitude and varying slowly with Z and Å is 
difficult to separate from other causes for dependence of Py on Z and A and from 
experimental uncertainties. From this fact and the assumption that the alpha particle 
formation probability does not depend especially strongly on the nuclear deformation, 
we obtain an explanation of the very reasonable values which one obtains for the 
nuclear radii from Po (Z, E4) by means of the simple Gamow theory. With the same 
assumption, the large factor (VIII-4) estimated by HimL and WHEELER?”” and amount- 
ing to — 40 (cf. Table 3) should have been easily detectable in the behaviour of the 
function Po (Z, EQ) and the nuclear radii calculated from it in the usual way. 
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b. Even parity alpha transitions in even-even nuclei. 


Most of the alpha groups which have been found in even-even nuclei are of 
the even parity type for which the nuclear transitions are 0+ — 0+,2+,4+4,.... 
The alpha groups corresponding to these transitions are denoted by Xp, %2,%4,..5 
respectively. The transition probability Ppo,1 for the group &; is given by the formula 
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which is obtained by specializing (VI-13) to the case that 1; = Ky = 0 (which implies 
that the alpha wave function on the nuclear surface contains only the angular moment- 
um component Q =— 0). Instead of using P%, 1 itself, we prefer to use the corresponding 
reduced transition probability c; which we define as the quotient of Po,7(E1) and 
P0,0(E1), where the explicitly denoted dependence on Ey; indicates that not only Pe, i 
but also Po,0 shall be calculated for the transition energy E; of the group &;. From 
the formula (VIII-10) we then get 
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if we note that it is only through the function (o (E7, Ro) that the expression (VIII-10) 
for Po,1 depends essentially on E;. It is useful to carry the reduction of the alpha 
transition probabilities for even-even nuclei one step further and to separate away 
not only the energy dependence, but also the /-dependence which is due to the ordinary 
centrifugal barrier. To this purpose, we introduce, instead of cz, a quantity b; which 
we define by 


… LARVER | 
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From this definition and the theoretical formula (VIII-11) for cz; we then get 
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For the quantities k;,v (B) corresponding to even values (0, 2, 4, 6, 8,...) of 
, ÆG . . . . . 
both I and I' we give the following numerical results, written in matrix form: 
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kK(0) = the unit matrix, 


ik] 
0.55 
Os 
0.01 
0.00 


k (1) = 


k (2) = 


i i: 
GS 01] 
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I = 2 
— 1.41 
1.51 
02 
0.47 
OG 


— 0.84 
1207 
OT 
OS? 
0:05 


— 0.41 
0.89 
=1D-5% 
0.06 
(ell 


=4 
0.64 
102 
1.30 
0595 
0.45 


0.28 
0601 
0.99 
— 0.60 
0.21 


0.07 
032 
0.89 
=— (31 
0.06 


res kr,v (0) = Or 
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0.59 
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From the empirical data on the alpha decays of even-even nuclei, PERLMAN and 
co-workers obtained the hindrance factors which are listed in Table 2. In this table, 
we have also listed the corresponding empirical c;-values which are the reciprocals 
of the hindrance factors, and the b;-values which we have obtained from the empirical 
cu-values according to (VIII-12), assuming Ry to be given by (VIII-7). The empirical 
knowledge thus obtained about the b;-values for the even parity alpha groups of 
even-even nuclei is summarized in Table 4. It is seen that b2 decreases as Z increases, 
ba has a sharp minimum for Z = 96, and b, has a broad maximum at about Z = 96. 
Also bg is rather large for this Z-value, but there is not much known about its vari- 
ation with Z. We shall now discuss this empirical evidence in connection with 


TABLE 4. Survey of the empirical b;-values for the 
even parity alpha groups of even-even nuclei. 


Parent 
nucleus ce by be bs 
(| 
Ra 1.2 —1.4 | 
Th 1.0 —1.5 0.7—0.9 04 
Wi 1.0 —1.3 0.6 —0.8 <0:4 
Pu 1.0 0.2 —0.3 03 0.2 
Cm 1.0 0.1 0.3 —0.4 0.3 
GE 0.7—0.8 0.2 | 03 
Fm 0.6 | 0.3 02 


the theoretical formula (VIII-13). From this formula it is seen that the b;-values for 
the even parity alpha groups depend only on the coefficients ag, a2,44,... and on 
the matrix elements kz,v (B) corresponding to even values of I and I'. Because of the 
appearance of the sign for the absolute value on the right-hand side of (VIII-13), it 
is also clear that, even if B is assumed to be known, one cannot determine ap, a2, 44,... 
uniquely from the empirical b;-values (for even /). Therefore it seems most reason- 
able to proceed by making simple assumptions about the shape of the nuclear surface 
and the alpha wave function on this surface, and to see to what extent one can 
understand the empirical data on the basis of such assumptions. According to the 
formulae in Chapter VI, it is not possible to separate the effect of the nuclear shape 
from that of the alpha wave function on the nuclear surface, y9. If wo is assumed not 
to change sign, it is also seen that this function can be replaced by a constant if at 
the same time the values of the actual deformation parameters $, are replaced by 
certain effective values. In the following, we therefore assume wo to be constant, but 
we must remember the effective nature of the nuclear deformation parameters 7. 
We shall now make some simple assumptions about the effective nuclear shape and 
discuss the b;-values thus obtained in relation to the empirical b;-values. 

We first consider the simple assumption that the shape of the nuclear surface 
is due to a pure quadrupole deformation and that the alpha wave function on this 
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surface is constant. The coefficients a; are then equalltorzeroktortl=70 Fan di from 
(VIIM-13) and the numerical values of kr,0 (B) for even / we get the b;-values which 
are listed in Table 5 for different values of B. The f2-values corresponding to these 
B-values have been obtained by means of the formula (VIII-6). From Q4-values 
obtained from Coulomb excitation data it follows that, for the nuclei around Th and U, 
the actual deformation parameter $2 has values of about 0.25. Empirically there is 
not much known about the variation of 2 with Z in the region of the very heavy 
elements, but for the moment we shall assume that $2 increases slightly with Z. Hence, 
we shall compare the empirical b;-values in Table 4 with those theoretical b;-values 
in Table 5 which correspond to f2-values lying in the neighbourhood of 0.25 and 


TABLE 5. Theoretical values of b2, ba, bø, and bg for the case that 
the nuclear surface possesses a pure quadrupole deformation 
and the alpha wave function is constant on this surface. 


P> | B bs by | be bg 
=():35 == 0.76 | 0.34 0.11 | 0.03 
= 0.4 2 0.61 | 0.20 0.05 0.01 
072, ER De ES SE ene MED 

0 0 0 0 | 0 | 0 

0.12 1 0.49 0.10 | 0.01 | 0.00 

0.24 | » 0.98 0.38 | 0.09 0.02 

0.35 3 1.35 0.74 0.27 0.08 


increasing slightly with Z. The calculated b2-values are of the right order of magnitude, 
but their slight variation with f2 is opposite to and less than what we should expect 
from the variation of the empirical b2-values with Z. The calculated ba-values have 
very roughly the right order of magnitude, but do not show the observed variation with 
Z. In fact, the calculated ba-values increase slightly with 2, whereas the empirical 
ba-values change fairly rapidly with Z and have a sharp maximum for Z = 96. The 
calculated b,-values increase slightly with 2, whereas the empirical b,-values have 
a broad maximum at about Z = 96. Furthermore, the calculated b,-values are smaller 
than the largest empirical b,-values. The bg8-values are empirically known only for two 
isotopes having the atomic numbers Z = 94 and Z = 96, and these bg-values are consi- 
derably larger than the small bg-values which we have calculated. Summarizing these 
results we see that, on the basis of the above simple assumption, we get the right order 
of magnitude for b2, whereas the higher b;-values deviate the more from the corre- 
sponding empirical values, the higher I becomes. Furthermore, the empirically found 
variations of b; with Z are not reproduced by the calculated b;-values for any l. 

We next assume the shape of the nuclear surface to be due to both a P2- 
deformation and a P4-deformation. The alpha wave function on the nuclear surface 
is still assumed to be constant. In a similar way as we have previously characterized 
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the P2-deformation by the quantity B, defined by (VI-9), we now characterize the 
P4-deformation by a quantity C which we define by 
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In the parenthesis on the right-hand side of this expression we may approximately 
neglect the second term. By using (VI-9) and (VIII-6) we then get the following 
approximate connection between C and fa: 
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According to (VI-10) and (VIII-14), we can express yp1 (Ø) in terms of C as follows: 
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For comparatively small values of C (e.g., for C £ 1), the b;-values can now be cal- 
culated by substituting the series expansions (VIII-16) into the theoretical formula 
(VIII-13). If 22 is assumed to be equal to the actual deformation parameter 82 which 
we, for the sake of simplicity, assume to have a constant value of about 0.24 (cor- 
responding to B = 2) throughout the whole region of the very heavy nucleif, and if 
Pa (= C/17) is determined such that the calculated b4-values agree with the empirical 
ba-values, we find the result shown in Table 6. Actually, there exists more than one 
possible C-value for every Z-value, but we have more or less arbitrarily required 
that (& shall not exceed 1, and, on this assumption, the determination of the C-values 
(in Table 6) from the empirical b4-values is unique. The magnitudes of the calculated 
bz-values as well as their decrease with increasing Z in the region from Th to Cm 
are in reasonable agreement with the empirical b2-values, but the slight increase with 
Z of the calculated b2-values in the region from Cm to Fm is not in quite satisfactory 
agreement with the behaviour of the empirical b2-values in this region. The calculated 
be-values differ from the empirical b,-values both as regards the orders of magnitude 
of these values and their variation with Z. It is also seen that the calculated bg8-values 


=£ For our later qualitative conclusions it is quite unessential whether /, is here considered to be 
constant or changing slightly with Z. 


NT. 3 es 45 


TABLE 6: 
Values of C, Ba (>= C/17), b2, be, and bg obtained by adjusting 
C such that the calculated b4-values agree with the correspond- 
ing empirical values. Here B has been assumed to be constant 
and equal to 2 (which corressponds to 2 = 0.24) for alithe nuclei 


considered. 

Parent by | 

nucleus (emp.) c Pa b» be bs 
ER 0.8 0.90 0.053 ets 0.39 0.16 
U 07 0.67 0.039 HLA 031 012 
Pu 0825 | — 0.28 — 0.016 0.92 0.02 0.00 
Cm OST | — 0.64 — 0.038 0.86 0.06 0.01 
Cf 02 — 0.40 =— 0.024 0.90 0.01 0:01 
Fm 0.3 — 0.16 — 0.009 0.95 0.05 0.00 


for the parent nuclei of Pu and Cm are appreciably too low. By considering both 
B2(= B/8.5) and p4 (>= C/17) as adjustable parameters, and determining them such 
that the calculated values of b» and b4 shall agree with the corresponding empirical 
values, we have found the figures in Table 7. The table shows that the effective de- 
formation parameter f2 is constant in the region of daughter nuclei between Ra and 
Pu, but decreases with increasing atomic number in the region of daughter nuclei 
above Pu. In connection with this result it may be noted that calculations of the nuclear 
quadrupole deformations indicate that the actual deformation parameter f2 increases 
slightly as one goes from Ac to Am, but may perhaps begin to decrease with increasing 
atomic number in the region above Am (private communication from Dr. NILSSON). 
It must, however, be borne in mind that our effective deformation parameter 2 may 
differ from the corresponding actual deformation parameter for the nuclear surface. 
The b46- and bg-values in Table 7 behave essentially in the same way as the corre- 
sponding values in Table 6, and they deviate considerably from the empirical values. 
Hence, it is not sufficient to assume only the two deformation parameters f2 and f4 


TABLE 7. Result of adjusting the parameters B(= 8.5 82) and C(= 174) such 
that the calculated values of b2 and b4 shall agree with the corresponding 
empirical values. 


De | 
EN ele GR 0 | R m | | 27] mn 
| | | | 

fie" PD) 0.26. | OS 0.041 URE SEERE 0578 ER 030 0.14 
U PD) 0.26 055 OOA RET DA 0 68 (ESP So tel 
Pu D19, 096 RER 04 == )/024 ØE | or 
Cm 235) 0.26 EO 0 047 RE 0:96 0.11 | 0.06 | 0.00 
cf 125 0.18 0.0 0.000 0.73 0.22 0103 007 
Fm 1% DE arr 0.024 0.65 0.30 0.09 | 0.02 
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to be different from zero. A rough estimate indicates that $,- and fØs-values of the 
order of magnitude of 0.02 should be employed to account for the large b,- and bg- 
values which have been found empirically for the parent nuclei of Pu and Cm. 

By means of Coulomb excitation experiments one could (at least in principle) 
obtain further information on the higher order multipole deformations which have 
been considered in this section, and also on the possible nuclear pear-shape which 
will be considered in the following section. 


Fig. 3. Two possibilities for the shape of a strongly deformed nucleus. 


Thexcircle"representstr-=rRQ 
The full-drawn curve represents r = Ro (1+0.25Y2 0 (9”)) 
The dotted curve represents r = RQ(1+0.25Y2,0 (9) +0.05 Y4,0 (9")) 


c. Odd parity alpha transitions in even-even nuclei. 


The even parity alpha groups %9, &2,... which we have discussed in the previous 
section are characteristic of all the strongly deformed even-even alpha emitters. For 
some isotopes of Th and U (and for Cm?2?) one has found also odd parity alpha 
groups 52, 53 which correspond to the nuclear transitions 0+—1—,3—,... and which 
will be denoted by ø1,%&3,..., respectively. The knowledge available on these odd 
parity alpha groups and on the corresponding odd parity states in the daughter nuclei 
is collected in the Tables 8 and 9%. The relatively small excitation energies of the 1— 


£Note added in proof: More recent data on the alpha groups populating the 1— states found in some 
even-even nuclei are given in the forthcoming review article on alpha decay by PERLMAN and RASMUSSEN232 
(see their fig. 15). In addition to improvements of the data in our Tables 8 and 9, they give data on the 
0+— 1— transitions of the parent nuclei Ra??? and RaZ?24, 
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TABLE 8. Data (from STEPHENS et al.) on the 1- states which have 
been found in some even-even daughter nuclei. 


Excitation energy | Hindrance factor 
É Rd of the 1- state of the &-group Degree of 
Nann CESRYDBR in the daughter going to the certainty of the 
parent nucleus = å 
nucleus 1- state in the 1- assignment 
(kev) daughter nucleus 
NL sko 242 i 4.0 certain 
Th?28 217 10 certain 
AN Sa 253 40 certain 
17% 232 ill certain 
ED SED 326 70 almost sure 
Cm2= 605 1000 somewhat doubtful 


states indicate that these states have the same individual particle configuration as the 
even parity states. There is also other evidence? that the odd parity alpha groups 
leave the daughter nucleus in a rotational band with Ky = 0 and that therefore the 
corresponding alpha particles are formed with 2 = 0. The odd parity states may be 
due to a soft asymmetric vibration of the type f3Y3,0(9") in the shape of the nuclear 
surface. (See, for instance, the review article on Coulomb excitation by ALDER et al.1.) 
Å special case of such a vibration is that of a nucleus which is pear-shaped in equili- 
brium and which may oscillate between the mirror shapes. Because of the rather 
long period for the vibrational motion it appears possible to consider the deformation 
parameter 3 as approximately constant during the penetration of the alpha particle 
through the interior part of the potential barrier. Since the shape of the nuclear sur- 
face is expected to be more important for the alpha decay than the anisotropy of the 
electrostatic field, the formula (VIII-10) can be used approximately also for the odd 


TABLE 9. Empirical b;-values for the odd parity 
alpha groups of even-even nuclei. 
(The b;-values are taken from Table 2). 


Parent 


nucleus a bs 
ll 

Th228 0.35 | 0.23 

Thz2e D7 | 

uæ 0.33 

U2 033 

Cm” 0.035 


parity alpha groups, even though it has been derived on the assumption that the 
differential equations for the alpha penetration problem do not couple the alpha 
groups going to different bands of the daughter nucleus. By defining c; and by for the 
odd parity alpha groups in the same way as for the even parity alpha groups we 
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then see that also the formulae (VIII-11), (VIII-12), and (VIII-13) apply to the odd 
parity alpha groups. The quantities k;, v(B) which appear in these formulae are dif- 
ferent from zero only if I and I' are either both even or both odd”. For odd values 
(SR SME) Kol bot bland V we give the following numerical values of kr, v (BB) 


written in matrix form: 


æn —3 [55 Ed 

0.61 = 068 0.37 014 

| 0.68 1.20 — 0.92 0.44 

EGE 0.37 = 0.92 1:28 = 0294 
£0) 74 0.44 — 0.94 128 

0.62 0:49 0.19 0.05 

— 0.49 0.97 — 0.59 0.21 

KEN" 40:19 — 0.59 0.99 — 0.59 
—0.05 0.21 — 0.59 1.00 

0.73 — 0.30 0.05 £ 0701 

| 0.88 LEN 0.06 
ken 005 —0.31 0.89 50131 
= 0701 0.06 50:31 0.89 


K(O)FEther una er (0) ES 07 


. . . . . . . . . . . . . . . . . . . . . . . . . . 5 
. . . . . . . . . e . . . . . . . . . . . . . . . . . . .: 
. . . . . . . . . ;= . . . . . . . . . . . . . . S a a = a 
. . . . . . . . . . s . . . . . . . . . . . . . . . . . . . 


1.60 0.58 0.10 0.01 

0.58 1851 0.53 0.10 

k(1)=4%. 0.10 0.53 1.48 0.52 
0.01 0.10 0.52 1.48 

DES IP 1.85 0.64 0217 

1.85 280 1570 0.58 

k(2) =3.. 0.64 1.70 DT 8 1.66 
OST 0.58 1.66 PETE 
(85:82 4.83 2839 0.79 

4.83 6.42 4.59 2.14 

k (3) = ) DA 4.59 6.02 4.39 
| DR) 2.14 4.39 5.92 

1=1 1=3 1=5 1=7 

£ Therefore the bj;-values for the even parity alpha groups depend on dyr do gone DUT TOON 


d1, d3, d5, + + »» Whereas the b;-values for the odd parity alpha groups depend on both these two sets of aj-Vvalues. 


Nr3 i 49 


By assuming the alpha wave function on the nuclear surface, yo (9), to be con- 
stant and by using (VI-10), (VI-11), and (VIII-13), we obtain to first approximation 


mm | ks) 2 | ns (BD fj] /Raf, —ERe) 
Kko,0 (B) ao ko,0 (B) å | x sd 
and 
bs k3,3 (B) 
Br ar DD] 


| 68] being understood as an average value of | 63 É Substituting appropriate numerical 
values in these two formulae, we get 


|63| > 0.051 

and 

bs 

55: 15% 
From the estimate of | B3| and the b3-values in Table 9 we find that | 63] - 0.01 to 0.03 
for the daughter nuclei corresponding to the parent nuclei of Th and U, and that | 65 | 
=> 0.002 for the daughter nucleus corresponding to the parent nucleus Cm?2?, The 
estimate of b3/b1 can at present be tested only for the parent nucleus Th??8 for which 
b3/b1 = 0.23/0.35 = 0.66 according to Table 9%. For other parent nuclei for which the 
alpha group &1 has been found, the estimate af b3/b1 may give information about the 
possibility of detecting the alpha group &«3. 


d. Favoured alpha transitions in odd-A nuclei. 


An especially interesting type of alpha decay is the favoured alpha decay which 
is Characterized by the alpha particle being formed from paired nucleons moving in 
orbitals differing only in the sign of the angular momentum component along the 
nuclear axis!5, In this type of alpha decay, the alpha particle is therefore formed with 
the angular momentum component 2 = 0, and hence the daughter nucleus is left in 
a band with Ky = K;. The alpha decay of an even-even nucleus is of the favoured 
type. The same is also true for the alpha decay of an odd-A nucleus if the last odd nu- 
cleon moves in the same orbital in the daughter nucleus as in the parent nucleus. All 
the favoured alpha groups of odd-AÅ nuclei which have been identified have even parity. 

Nucleons occupying states which differ only in the sign of the angular momentum 
component along the nuclear axis interact especially strongly due to the large overlap 
of their wave functions. Hence it is most probable for an alpha particle to be formed 
from such paired nucleons, i. €., favoured alpha decay is intrinsically more probable 
than other types of alpha decay. The actual alpha transition probabilities are, however, 


£ Note added in proof: Quite recent, unpublished data on odd parity alpha transitions in even-even 
uclei indicate similar discrepancies between the theoretical and empirical values of b,/b, to those for the 
arent nucleus Th??8 (private communication from Dr. MOTTELSON). 
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strongly affected by the potential barrier, and in many cases they are therefore much 
smaller for the favoured than for some unfavoured alpha groups. 

Because of the large intrinsic probability for favoured alpha decay, the favoured 
alpha groups of even parity emerging from an odd-A nucleus are recognized by the 
fact that their F-values (see Chapter VII) are comparatively large, being less than 1, 
but still of this order of magnitude for small spin changes. Thus, the F-value for the 
special alpha group of this type which corresponds to the nuclear transition involving 
no change of spin is found to be about 0.5 on the average (see Table 1). We may 
express this fact simply by saying that for odd-A nuclei the favoured alpha groups 
of even parity are intrinsically hindered by a factor of about 2 compared with the 
corresponding alpha groups for even-even nuclei. 

For favoured alpha decay the alpha wave function on the nuclear surface, 
wo (9”), contains only the angular momentum component Ø— 0 and is therefore 
axially symmetric with respect to the nuclear axis. This fact implies that the coef- 
ficients a, 9 defined by (VI-11) are different from zero only for Q = 0. For favoured 
alpha groups the formula (VI-13) therefore becomes 


[ Ro 
> ev 
I Kf,14 ng Go (Er, ka) E 


Furl MeD Ea | 
0 


SKE EET ER 
(VIII-17) 


We define the reduced transition probability cxr,74 for the favoured alpha 
transition (K;, I;) > (Ky, I4) as the quotient of the probability for this favoured alpha 
transition and the probability for the favoured alpha transition which involves no 
change of spin and parity, where the latter probability is corrected such that it corre- 
sponds to an alpha transition energy equal to that for the transition (K;, I;) — (Ky, I4). 
From this definition and (VIII-17), and the fact that the essential energy dependence 
in this expression for Pky,14 enters through Go (E14, Ro), we get 


2» ERR OT IEL BSK) 
1 


CKy, ly (VIII-18) 


Dre SER OPRET KDE 
I 
where 


Jexp)- HD NE 4 ii (VIIT-19) 


In the numerator of (VIII-18) the sum over I is restricted either to only even /-values 
or to only odd /-values depending on whether the parities of the parent and daughter 
nuclei are the same or the opposite (cf. the last few lines in this section). In the deno- 
minator of (VIII-18) the sum over I is always restricted to only even /[-values. 
The favoured alpha decay of odd-A nuclei must in many respects be very 
similar to that of even-even nuclei. It would therefore be natural to assume that, 
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for favoured alpha decay of an odd-A nucleus, the alpha wave function on the nuclear 
surface, wo (9”), could be calculated by interpolating the corresponding functions 
yo (9”) for neighbouring even-even nuclei to the Z-value of the odd-A nucleus con- 
sidered. However, according to the formulae (VIII-10) and (VIII-17), this assumption 


TABLE 10. Data on those favoured alpha decays for which one can test the formula 
which relates the relative intensities of the favoured alpha groups of an 
odd-A nucleus to the relative intensities of the alpha groups of neighbour- 
ing even-even nuclei. 


Parent Fe: es Er CK+,I 4 CK£ I 
KF=K=1L I eg FSK i KEN A 
nucleus I (MeV) emp. ERNE REM je theor. 
STORA 5727 1 1 
91Pa2?31 3/2 SP VA 0.32 i 0.87 | 0.091 | 0.38 
7/21 /84:593"(cale) | 0.23 
| 5/2 | 4.8157 il 1 
U?2334+ 5/2 17 0 047732 0.34 1 0.83 | 0.071 0.31 
| SDR NEL 7474 0.093 0.13 
| 11/2 | 4.650 (calc.) 0.018 
|| 
oOo kel PM EA Før Pa 
Ko 550 1 1 
94Pu?3? 1972 | 39: 15137 0.33 1 0.65 | 0.011 0.0026 0.26 
(5 ro RES 099 0.35 0.39 
7/2 | 5.070 (calc.) 0.0049 
ma Re ars Fa 
(85054820 1 i 
HE 7 07 15 4391 0.31 0.24 
ssAm241 5/2 | 9/2 | 5.3860 0.083 il 0.59 | 0.0040 | 0.0031 0.083 
(ry 55270 0.0018 0.0016 
DE orn 0.00076 0.0013 
cz N KHÆSAETS — 
Bo 267 1 1 
gå mis 5/2 TJ MED DDA 0.25 1 0.59 | 0.0040 | 0.0031 0.24 
9/2 | 5.169 0.052 | 0.083 
| 11/2 | 5.102 (calc.) 0.0016 
ss fENnE zd : my 
[52 51777 1 i! 
gm 5/2 HG JR 55732 0.29 il 0:570100:00177 700033 0.23 
9/2 5.674 (calc.) 0.080 
7/8 CV" 67633 (beg | i 
(ed DE 72 9/2 | 6.592 (0:18 1827 50:290 0 01400022 OMI 
| 11/2 | 6.545 0.050 | 0.033 
| 13/2 0.0034 


Note added in proof: Recently, GoLnin et al." have investigated the decay af U?%? and observed two 
new alpha groups corresponding to transitions to the states of the daughter nucleus with spins 11/2 and 
13/2. The intensity values, which are given with large uncertainties, are only in qualitative agreement with 


the formula (VIII-18). 


7% 
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does not account for the previously mentioned hindrance by a factor of about 2 for 
the favoured alpha groups of even parity in odd-A nuclei. Therefore we make the 
somewhat more complicated assumption that, for favoured alpha decay of an odd-Å 
nucleus, the function wo (9”) is equal to the function which one obtains by first inter- 
polating the corresponding functions yo (9”) for neighbouring even-even nuclei to the 
odd-A nucleus considered and then multiplying the resulting function by an angular- 
independent factor having an absolute value of about 0.7 (= V0.5). According to 
(VI-10), (VI-11), (VIII-11), and (VIII-19), the quantities c, for our odd-ÅA nucleus 
are then proportional to the reduced alpha transition probabilities c; for even-even 
nuclei, interpolated to the odd-A nucleus considered. In the formula (VIII-18) we 
may therefore simply replace the c;-values by these interpolated reduced transition 
probabilities cz. In this way, we get an intensity formula which makes it possible to 
calculate the reduced transition probabilities for the favoured alpha groups of odd-A 
nuclei on the basis of the empirically determined reduced transition probabilities c; 
for the alpha groups of even-even nuclei. In a slightly different form, this formula 
was first given in a paper by BoxHr, FRromaAnN, and MOTTELSON?5, There, this intensity 
formula was tested by comparing, for each possible odd-A parent nucleus, the relative 
values of the empirical transition probabilities for the favoured alpha groups with 
the relative values of the corresponding theoretical transition probabilities. A more 
sensitive way of testing the formula is, however, to compare the empirical values 
of cx,,1p4 for the favoured alpha groups of the odd-A nuclei (obtained from Table 1) 
with the corresponding theoretical values. This comparison is given in Table 10 which 
lists those favoured alpha groups of odd-A nuclei for which sufficient empirical data 
are available. The values of Kr (= K;) and I; are taken from the Tables III and IV 
in the previously mentioned paper!5, except for the parent nucleus E253 for which 
the suggested identification of the favoured alpha groups is new. For this nucleus 
we have assumed that K/; = 7/2. The values of the alpha particle energies E,, and the 
empirical values of the reduced transition probabilities CKp,Ip are taken from Table 1 
in the present paper. The cz-values are obtained by interpolating between the empirical 
crvalues for the even-even nuclei (see Table 2). The theoretical values of cx/,174 have 
been calculated from (VIII-18) by neglecting the terms corresponding to 1>6 (and 
in most cases also the term corresponding to I = 6). When the empirical and the 
theoretical values of cxyp,174 in Table 10 are compared, one finds an agreement which 
is as good as can be expected considering the approximations involved in the theore- 
tical formula and the uncertainties in the experimental data. 

Although no favoured alpha groups of odd parity have so far been identified 
for odd-A nuclei, the existence of odd parity alpha groups for some of the heavy 
even-even nuclei suggests the possibility for the existence of favoured alpha groups 
of odd parity for neighbouring odd-A nuclei. In the case of such an odd-A nucleus, 


x Note added in proof: The identification of the favoured alpha groups of E253 js also given in the 
forthcoming review article on alpha decay by PERLMAN and RASMUSSEN?" (their table XVI). Furthermore, 
these authors suggest an identification for the favoured alpha groups of Pu2?t, 
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the two bands in which the favoured alpha groups of even and odd parity, respectively, 
leave the daughter nucleus are expected to have similar level structure, but opposite 
parity, and to be displaced with respect to each other by roughly the same amount 
as the even and odd parity bands in neighbouring even-even nuclei, i. e., by a few 
hundred keV. The possibility for detecting the odd parity alpha groups here suggested 
can be estimated on the basis of the intensity formula (VIII-18) and the empirical 
c1-values for the even-even nuclei. 


e. Angular distributions of the alpha groups from polarized nuclei. 


The expression (VI-14) for the angular distributions of the alpha groups emerging 
from polarized nuclei can be evaluated by means of the formula for the expansion 
of a product of two spherical harmonics in terms of spherical harmonics and by 
means of a sum rule!!a which expresses sums (over magnetic quantum numbers) of pro- 
ducts of three Clebsch-Gordan coefficients in terms of Racah coefficients. Thus, if the 
coefficients aa; p åre assumed to be real, the expression (VI-14) becomes 


SÅ ra REE brE> ne LT De BERETTE) | 
Å i; er (VIII-20) 


x SC (L=-TÆL) cCoSs KEDE GE Se De ARIG TET RENEE DER EN | 
KER 


7 
where 
2 1(1+1) 
ig = expl- ( a - vær 


SED? (GK 0-02] I, LI, KODER) org 
- 


| (VIII-21) 


and the coefficients Z are related to the Racah coefficients W according to the formula!!a 
ZA (USE DEACON | d) 
=(- i rbee2 AE a+ 1 Q2b HJE DØDT? INN Gordo ge es gf Mage | (RER OD)E 


The expression (VIII-20) is to be averaged with respect to the initial distribution of 
M;-values. This averaging affects only the Clebsch-Gordan coefficient which appears 
in (VIII-20). For L = 0, this Clebsch-Gordan coefficient is always equal to 1, and for 
L = 2 its average value with respect to the initial distribution of M;-values is 
MER USD ANSE Dig Al. AVE ZN 
IkSDED Pa ØE fy MIS LAVE = == ere å i =22 
GØ SR RE Por EDEL] 
where Å is defined by 
1 
A SE ME —41;(1;+1) j (VIII-23) 
æl G i AV 
i 
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and is a measure of the polarization of the parent nuclei, being equal to zero for non- 
polarized parent nuclei. Considering only even L-values and using (VIII-22), we can 
write the above mentioned average value of the expression (VIII-20) as follows: 


2 
Ro DYBE 
(Å F; 
ASA (ens, 1 2 Vine 4 
VIF FR ERKED ng i BETA ) 2( ) f ( ) 


where n 
(DES VS : 
VETED ROED OT ORE) 


C14, 14 (B) HÆR 


ES SE DRS se DS ETERN ANSER NEN 3] 
BENS 
l 


The case of favoured alpha decay of parent nuclei of spin 5/2 is especially 
interesting for the interpretation of the experimental results available on polarized 
nuclej?2, 46, 46a, 47, Assuming the function y1 (9”) to be independent of 9” (i. e., all the 
coefficients 7, except ap, to be zero), we obtain the numerical values of the coefficient 
C. s listed in Table 11. The quadrupole deformations found for the very heavy nuclei 


2, 


tølg 


TABLE 11. 
Values of C6 5 (B) for the case of favoured alpha decay 
corresponding to x = 50, kRo/x = 0.2, and y1 (9”) being inde- 
pendent of 9". 


É C3,5 (B) 
=—2 26 
el 0:83 
0 0 
1 1.15 
2 2.07 


correspond to prolate nuclear shapes, i. e., to positive B-values. For such B-values, 
the coefficient C; , is positive, while the experimental results appear to indicate an 
alpha distribution peaked perpendicular to the direction of polarization?2, 46, 46a, 47, 
We abstain from giving values for C6 2: 
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Appendix Å. 
Formulae and Graphs for the Gamow Wave Functions. 


In this appendix, we collect some formulaef for Gamow's well-known WKB- 
solution of the differential equation describing the penetration of an alpha particle 
of given angular momentum I" through a spherically symmetric Coulomb barrier. For 
1=0, we also give graphical illustrations of some properties of this solution. 

Consider the differential equation 


2 2 


KDE EK ED) 


ie 2 mr? 


ET (TER 


where 


UKE) 


The Gamow wave function G;(E,r) is a solution of (A-1) representing an outgoing 
wave for large values of r. Such a solution is unique except for an arbitrary constant 
factor, and by choosing this factor conveniently, we get 


æn Nes 
Er 4 i . es, | 
| E ) esplæf Sr; | Di when U;>E 


ha gig 6 . Rg: 3 
E 5 esp lig +ikl (2 dt) ar] when U;<E 


Gi (E, T) ==; (A-2) 


TE | E 


according to the WKB-approximation, which in general is considered to be valid if 
the condition 


<<1 (A-3) 


is fulfilled. The two integration limits in (A-2), which are not written out, both corre- 
spond to the point where U;, which is an abbreviation of U;(r), is equal to E. The 
quantity k which appears in (A-2) and (A-3) is defined by 


1 -— 
k==/2mEeE, 
mL 
and represents the magnitude of the wave vector for large values of r. The corresponding 
væ E 
v= | . 
m 


For the evaluation of (A-2) it is convenient to introduce the dimensionless 


velocity is 


quantity 


+ Many of these formulae are already known and can be found, e. g., in the book by GAMOW and 
CRITCHFIELD28, 
8+ 
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AGE) E 
hiv 


= (A-4) 
which is well-known in atomic collision problems. In our case, x is of the order of 
magnitude of 50. The fact that x>>1 is a criterion that a classical orbital picture 
is possible outside the potential barrier!”, 

For the quantity (U;— E)/E which appears in (A-2) and (A-3), we now easily 
find that 


5 i 1 
ve (2 en (A-5) 


E |kr 


(0) 


(9 05 1.0 15meker 
Fad 


Fig. 4. Graphical representation of the expression on the left-hand side in the condition (A-6) for the validity 
of the WKB-approximation. 


In the condition (A-3) it is for our purpose a sufficiently good approximation 
to replace U; (r) by UQ (1). According to (A-5), this is the same as replacing (U; — E)/E 
by x/kr—1. In this way the condition (A-3) becomes 
k kr |3Y2 1 
(kr 8% ) SEE (A-6) 
i x % 
A graphical representation of the expression on the left-hand side of this inequality 
is shown in Fig. 4, from which it is seen that the condition (A-6) is fulfilled except 
for a small region around the classical turning point. 
In the following we shall treat the two cases kr<x and kr>x separately. 


i 
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Instead of r we introduce the variable & which we define by 


eN E 'z værer A-7 
x U4 Gt) OSTE (4 > a) 


Es 


We note the two formulae 
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and 
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for"later" use: 
The condition (A-6) for the validity of the WKB-approximation becomes 


; 1 
SIN COSIK=——% (A-10a) 
x 


From (A-2) (for U;7>E), (A-5), (A-7a), and (A-10a) we get the well-known ap- 
proximate formula 


GE) E= (cot x)" exp É («-sin a cosa) + Dig af (A-11a) 


which is valid if the centrifugal barrier is small compared to the Coulomb barrier. 
We shall now calculate G;(E+4E, r)/G1 (E, r). Obviously, 


G(E+AE, r) mic i SR | 
En LE SEE | 
Using this formula, together with (A-8a), (A-9a), and (A-11a), we get 


re as TT å 1(1+1) 
GÅEz) SN 1 x(x+ sin&cos&)+ SEG El 


if the condition (A-10a) is also taken into account. 
The formulae (A-11a) and (A-12a) are valid only if it is assumed that the centri- 
fugal potential is small compared to the Coulomb potential. If this condition is 


In AB: 


not introduced, one finds, by a straight-forward evaluation of the integral appearing 
in the WKB-approximation formula (A-2) for U;>E, that?3 


1 


» De Roe ae bra ET NR 
G e | -A årrslves xp (SE A-13: 
BCBST) 2 i) Fr 1! exp Æ LES bs i ( a) 
where 
x+2a+ 1 +æ — ax? WESER 
F(a,x) =|/aln RA, V4a(a+ mær. Lar COS SE SER ERE HOS LERES) 
cV1+4a Vare 


and arccos shall be chosen to have a value between 0 and ax. If we suppose I (1+1)/x2 
to be small and expand F(1(1+1)/x2, kr/x) in terms of powers of 1 (1+1)/x2, keeping 
only terms that are linear in /(1+1)/x?, we find from (A-13a) and (A-14a) the ap- 
proximate formula 


[(1+1 
Gi (ER) = G(E,r) esp | S GE æ (A-15a) 


which gives the /-dependence already found in (A-11a). 
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TABLE 12. For the case x = 50 and kRy = 0.2 x, values of the relative penetration 
probabilities for different angular momenta / are listed. They are given in 
two different degrees of approximation. The last column corresponds to 
column 3 in Table 13. 


iG, (E, Ry)/G1 (E, Ry))" iG, (E, Ry)/G1 (E, Ro))” 
l according to according to l according to according to 
(A-13a) and (A-14a) (A-15a) (A-13a)and (A-14a) (A-15a) 
0 1 1 8 0.386-10=2 OST5- 1052 
1 0.854 0.852 9 0.995-1072 OASE 
2 0.624 | 0.619 10 0.224:1073 OSTE RO re 
3 0.390 0.383 11 0.440:1077 0.259:1057 
4 0.209 0.202 12 05702105 0.380"1055 
5 ODO 05 OFPOFEO ET 13 OML6-105 0:475-1058 
6 OF 10 0.347:1077 14 WIE SETDES 0506510 
7 O30 0 IR OSS SKO al 15 OM90 FO 0.459:105£ 


Using formulae equivalent to (A-13a) and (A-14a), DEvanEY?3 has calculated 
the relative penetration probabilities (Go (CERRO ENER RO for alpha particles of 
different angular momenta / in the case that Z—2 = 86, Ro = 9.87:10713 cm, and 
FRA 88 MeV (er 4 OM Rand ER 0105) ESS numericalkresultsll(tablerlkon 
p. 591) which are also listed in BLaTT and WEI1sskorr?? (table 3.1 on p. 577) appear, 
however, to be considerably too low for large values of l. It is probable that DEVANEY 
has calculated (Go (E, Ro)/G1 (E, Ro) instead of (Go (E, Ro)/G1 (E, Ro)P”. For the 
parameter values x = 50 and kR/x = 0.2, which differ only slightly from those 


TABLE 13. Relative penetration probabilities for different angular momenta /, calcula- 
ted according to the simple approximate expression expt-cl(1+1)). The 
values chosen for the parameter c, which is related to the atomic parameters 
according to the formula (A-16a), cover the range of practical interest. 


exp (—cl (I + 1)) 


for e0:07 for 0140/08 for c = 0.09 | force 0M0 
0 1 1 jl al 
1 0.870 | 0.852 0.835 0.819 
2) 0.657 0.619 0.583 0.549 
3 0.432 0.383 0.340 0.301 
4 0.247 0.202 0.165 0135 
5 023 | 0.907:10—71 O6 72: KO 0.498:10—1 
6 0.528:1071 0.347:1071 09228105 0.150:1071 
EL 0.198:1077 OM sT 0 01647105 0:370:10- 
8 0.647-10—2 VS175-1052 OSSE 07475105 
9 01837105 0747-1052 0:304:708 OSPF OEE 
10 0.453:1073 0351-7052 0.502.107 0.167:10— 
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corresponding to the case treated numerically by DEVANEY, we give in Table 12 the 

values of Go (E, Ro)/G1(E, Ro)Y according to (A-13a)together with (A-14a), and (A-15 a). 
According to the simple approximate formula (A-15a) the relative penetration 

probabilities for different angular momenta / are expf(— cl (1+1)), where 


2 Ø m Dyr JONES DYE 

> 1 SEE. AS bra 
i E| kRo ÅL ER, ( ) 
and Ry is the radius of the spherical nucleus. In Table 13 numerical values of 
exp(— cl (1 DD are listed for the range of c-values which is of interest for experi- 
mental interpretations. 


Fig. 5. Graphical representation of the 10-logarithm of |G, (E, r)|. 


60 Nr. 3 


b. The case kr>x, i.e., Up (r)<E. 
Instead of r we introduce the variable 4, which we define by 


kr E 


cosh? A =— = (4>0). (AS4D) 


xx UQ(r) 


logIG (E+AE ri log IG(E,r)l 
A 


Fig. 6. The increase of the 10-logarithm of |G, (E, T)|, as E is 
increased by AE but r is fixed, is obtained by multiply- 
ing the appropriate ordinate in this graph by AE/E. 


We note the two formulae 
oa sh A 
£ DE Es (A-8b) 


and 


OA. coth AX 


GE BE CB) 


to be used later. 
The condition (A-6) for the validity of the WKB-approximation becomes 


2sinh? AcoshA >>. (A-10b) 
x% 


From (A-2) (for U;<E), (A-5), and (A-7b) we get the approximate formula 
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f1(1+1) 1 | 


Gr (EST) = (Cothla) PI 4x2  sinh2A cosh2A | (A-11b) 


+1 


T +x(sinh A cosh 4 — 4) - Aus: Liga) |. 


Now, we calculate | G(E+4E,r)/Gr(E,r) | . Obviously, 


T2 G(E+AE,r) 
GE 2) 


9 
DE In | GE, 7) 4£. 


Using this formula, together with (A-8b), (A-9b), and (A-11b), we get 


SGÅESA Er) f 1 LOESD NE LER 4] 
2 REED ER | te ik, EA | 49 
G(E,r) ØR (| ED I E-AR FAN 52545) 
log | Go(E + A ER Tr) Go (E-7)| 


From the formulae in this appendix, log | Go(E,r) | and - F TEIL 


can easily be found as functions of kr/x and x. Figs. 5 and 6 give graphical repre- 


sentations of these functions. In connection with these figures, it should be remarked 
that, in the approximation considered here, Go(E,r) is real and positive for kr<x, 
whereaskitkiskcomplexE for kkr ÆT heRSameRiskobviouslyktruerforkthertunceton 
G(E+AE, T)/G( (E, Tr). ; 


Appendix B. 
Penetration of an Alpha Particle through a Static Non-central Field. 


The wave equation for an alpha particle penetrating a three-dimensional potential 
barrier, which is" fixed in space, is 
aps KELans (B-1) 
where 


- 


See vær 
KK (re // gl) E)= k]/ = 


V denoting the potential, m the mass of the alpha particle, E its kinetic energy after 
it has escaped completely from the potential barrier, and k the magnitude of the 
corresponding wave vector. For the construction of an approximate solution of (B-1), 
CHRIstY has recently developed a three-dimensional WKB-method”. By means of 
such a solution, he can approximately express the alpha wave function in an arbi- 
trary point P inside the potential barrier in terms of the alpha wave function on 
the nuclear surface, except for a constant factor. In this appendix, we combine 
CxariIstyY's ideas with a method similar to that used by KiRCHHOFF for the rigorous 


+ Private communication from Professor R. F. CHRIstY to Professor A. BoHRr. See also reference 21. 
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justification of Huyghens” principle in optics" (see BorN's monograph on optics"8,) 
The result of these considerations is an approximate integral representation of the 
alpha wave function inside the potential barrier in terms of the alpha wave function 
on the nuclear surface, which is closely related to the result of CHristy although 
more rigorously justified. The kernel H(P, P') appearing in our integral represen- 
tation corresponds, in KirCHHOFrF's formula for light waves in free space, to the spherical 
wave explik| Tp TP 13/1] Tp -—rp]. Considered as a function of the point P", the kernel 
H(P, P') is a solution of (B-1) with a singularity behaving as 1/|rp-r, | ASTI 
We first consider the construction of this kernel by means of the second order ap- 
proximation of the three-dimensional WKB-method, then we derive the previously 
mentioned integral representation for the alpha wave function. To investigate the con- 
sistency of the resulting formula we apply it to the case of a spherically symmetric 
potential barrier and find that the well-known WKB-approximation solution of this 
problem is very nearly reproduced. Finally, we use our integral representation to 
obtain a simple approximate solution of the problem concerning the penetration of 
an alpha particle through an anisotropic potential barrier. 
By substituting 

== GP (B-2) 

into (B-1), we get 
[grad S FSAS= "KE, 


In analogy to the procedure in the one-dimensional WKB-method, we here neglect 
the term AS and get 


ferades re (B-3) 


This equation is similar to the eikonal equation in geometrical optics (with S/k 
corresponding to the eikonal) or to the Hamilton-Jacobi equation in classical mecha- 
nics (with AS corresponding to HAMILTON's characteristic function). From (B-3) we get 


gradS = K(r)es, (B-4) 


where es is a unit vector field; its determination will be discussed later. For an arbi- 
trary choice of the vector field es there exists in general no solution $ of (B-4). The 
necessary and sufficient condition for the existence of such a solution is that 
curl (Kes) = 0. Itis easily seen that, if there exists a solution S of (B-4), the line integral 
VK ds, taken along a curve which is everywhere tangent to es and which terminates 
on any two surfaces of the type S = constant, has the same value for all such curves 
between these two surfaces, and it is also seen that, for any two points P and P' on 


any such curve, this -curve is the path between P and P' which minimizes the 


ap" 
integral IV: ds. 
-P. 


An extension of Huyghens” principle in optics to the case that the index of refraction is not con- 
stant DE varies slowly in space can also be obtained along similar lines as those which are described in this 
appendix. 
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To construct a special solution of (B-3), we determine the path which joins two 
points P and P' and which has the extremal property 


AP' 

|x ds = minimum, (B-5) 
OG så 

s being the arc length measured in the direction from P to P". The unit vector which 
is tangent to this curve and which points in the direction of increasing arc length s 
is denoted by es. The requirement (B-5) is analogous to Fermat's principle in geome- 
trical optics, and the path which one obtains from it is analogous to a light ray. 
Therefore this path will be called the ray between the points P and P'. From the 
general formula 


aP' P' ig d | 
bå Kds= Key dr SÅ grad K—— (Kes); : år ds (B-6) 
ÉuP ze Ad ds j 


it follows that the differential equation for a ray is 
P d ”g 
grad K— E (Kes)=0, (B-7) 


for the first term in (B-6) vanishes for fixed P and P". By keeping P fixed, but letting 
P" take all possible positions, we can cover the space by rays which emanate from 
P. From (B-6) and (B-7) it is then clear that all these rays are cut orthogonally by 
any surface such that the points P" on it fulfill 


op 

| ds =-constant: 

vP 

where the integration is performed along the ray from P to P". From a simple geome- 
trical consideration it then follows that the condition curl (Kes) = 0 is fulfilled. We 
now define 


P 
S= N Kds, (B-8) 
wP 
where the path of integration is the ray between the given point P and a variable 
point P". This function $ is obviously a solution of (B-4) and therefore also of (B-3), 
if es on the right-hand side of (B-4) is defined as above by means of the paths obtained 
from the extremal requirement (B-5). Hence e75, with $S defined by (B-8), is a first 
order WKB-approximation solution of the Schrodinger equation (B-1). For the validity 
of this approximate solution it is sufficient that simultaneously 


1 
grad K SM] 
dives <1 
K iz . 
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The first of these two inequalities is identical to the well-known condition for the 
validity of the first order WKB-method in one dimension, whereas the second one, 
which, in the case of a spherically symmetric potential barrier is certainly fulfilled 
if r K (r) is large compared to 1 and if furthermore the point P' does not lie too close 
to P, appears only for the WKB-method in three dimensions. 

Next, we derive the second approximation of the solution of (B-1) which we 
have found above in first approximation. To this purpose, we write 


DE GT, CE) 
S being still defined according to (B-8). By substituting (B-9) into (B-1), and using 
(B-4) or (B-8), we get - 
u div (Kes)+2 Kes gradu—- Au=0. 
In analogy to the procedure in the one-dimensional WKB-method we neglect the 
term Au and get 


dl k 
es-grad In (uk?) rr Ndivesrt 


fø £ 


u= const ke exp | | idives dl 


the path of integration here being along a ray emerging from the point P. From this 
expression for u, and from (B-8) and (B-9), we find that the function 


PINE ap" rr ap D 
NEDE ED MEN DE 
HK DABA) le) exp ) É: 5 dives ST V2 (P') exp , | SA) (B-10) 


is an approximate solution of (B-1). Here, the point Pis considered to be fixed and 
the point P' to be variable, and the integrations are performed along the ray joining P 
and P" (see Fig. 7); P1 is an arbitrary point on this ray (conveniently chosen to lie 
close to P). Because all the rays go through the same point P, a simple geometrical 
consideration shows that 


1 


SÆL SEP: 
[re -rp]| as — P 


I dives — 
at the point 
tthetponters 


If P" lies sufficiently close to P (and if P1 is chosen to lie close to P), the formula 
(B-10) therefore becomes 


1 RK 
FHLGESBA) FF EEr 2] j KKP) exp LN (r) AN 
BE 
Hence, 
F: 1 
ITÆPI PS) > SE ASSER 


[ry —rp | 
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For the validity of the approximate solution (B—10) it is sufficient that simultaneously 


læst 1 

KkÅK FE grad , eæl 
er die; dive, 
|jeørad k K B]R 


The first series of these inequalities is the same as the condition for the validity of 
the second order approximation of the one-dimensional WKB-method, whereas the 
second series of the above inequalities appears only in the three-dimensional WKB- 
method. Ås P' approaches P the inequality | div es|< <K becomes much violated ; 
nevertheles, (B-10) may still remain an approximate solution of (B-1) and, hence, the 


Fig. 7. The figure shows rays emerging from the point P. On the ray PP", we have indicated the point P; 
and the unit vector es along the ray in an arbitrary point on it. 


singularity which we have found for H (P, P') may be relevant. Thus, a consideration 
of a spherically symmetric potential barrier indicates that for this case the approximate 
grad ele 
The function H(P, P') which we have constructed and discussed above will appear 
as the kernel in the approximate integral representation of the alpha wave function 
which we now derive. ! 

Considered as a function of P" the function H(P, P"), defined by (B-10), is a 
second order WKB-approximation solution of (B-1) having the singularity 1//rr—r 2] 
at 7p = rp. The alpha wave function w (P") inside the potential barrier is also a 
solution of. (B-1), and therefore we get” 

p(P) =—— NÅ fy (P') grade H(P, P')-H(P, P') grade p(P")) - dap (B-11) 


Å 77 


solution (B-10) is valid in the immediate neighbourhood of the point P if 


from GREEN's formula, the surface element at P" being denoted by døp'. In the formula 
(B-11), which is analogous to KiRkRCHHOFF's formula in optics, we consider the inte- 
+ The formula (B-11) is exact as far as the kernel H (P, P") with the required singularity at 7,7 = rr, 
is a correct solution of (B-1). As has been pointed out by GROENEWOLD?, it is possible in principle to use 
an exact formula for this kernel. 
Mat. Fys. Skr. Dan. Vid. Selsk. 1, no.3. 9 
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gration to be performed over the nuclear surface and over the surface of = 
sphere concentric with the nuclear surface and conveniently chosen somewhere in 
the potential barrier (see Fig. 8). Since the alpha wave function wp on this outer 
surface of integration is small compared to the same wave function on the nuclear 


z 

surface, and since H (P, P') contains the factor exp E K ds , the integral over the 
P 

outer surface is negligible compared to the integral over the nuclear surface, unless 


the point P lies in the vicinity of the outer surface. In order to obtain the alpha wave 


de 


Fig. 8. Geometrical illustration of the quantities appearing in (B-11) and (B-14). The integrals in these 
formulae shall be performed over the nuclear surface and over the larger spherical surface which lies 
somewhere in the potential barrier. The surface elements do” on these two surfaces have the direc- 
tions shown in the figure. Unless the point P lies in the vicinity of the outer spherical surface, the in- 
tegral over this surface is negligible compared to that over the nuclear surface. The curve PP” is the 
ray joining the point P in the potential barrier and a point P" on the nuclear surface, and es is the 
unit vector along this ray at the point P". The gradient of the alpha wave function at the point P" 
has the direction -æ where e is a unit vector which we assume to be approximately orthogonal to 
the nuclear surface at the point P". For the most important rays PP", the unit vectors es and e 
have approximately opposite directions. 


function in the interior region of the potential barrier, where we wish to apply this 
approximation, it is therefore sufficient to perform the integral in (B-11) only over 
the nuclear surface. To simplify (B-11) further we note that, in first approximation, 


STA dl PAT (ESRB O'-=ESKECRAPEN (LÆRER (B-12) 
according to (B-10). Furthermore, if p varies slowly on the nuclear surface compared 
with its rapid variation in the radial direction, we can write approximately 

grad pay (DP ) FE KELD (PI ler (B-13) 
where e denotes a unit vector which is directed almost radially outwards from the 
nuclear surface. The minus sign in (B-13) is due to the fact that p is decaying 


exponentially in the outward direction. Substituting (B-12) and (B-13) into (B-11). 
we get the approximate formula 
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p(P)= SÅ | H(P,P)K(P') p(P") (e, —e): dop', (B-14) 
where the integration is to be performed over the nuclear surface. The vector e is 
directed essentially radially outwards from the nuclear surface, whereas es and dop' 
are directed essentially in the opposite direction (cf. Fig. 8). The formula (B-14), 
which is analogous to a formula appearing in KircHHOFF's treatment of the diffraction 
of light"S, has a similar relation to CHRiIsTty's approximation method as the just men- 
tioned formula in KircHHorrF's theory has to HUuYGHENS' principle. 

To illustrate the use of (B-14) and to investigate the accuracy of this formula, 
we now consider the well-known case that the potential V(r) is spherically sym- 
metric. Using polar coordinates (r, å,9), we find that the extremal problem (B-5) 
can be written 


An ( 


(Pr p , NØ) 3 
Å K (7) Eg + r? rå + sin?d (SP) | dr = minimum. 
Ll Beg Å 


The corresponding Euler-Lagrange equations give 


K (7) r? sin? de 
E =e (B-15) 
fle 12 ; FE 2 4 = Re 
REE (7) + sin?4 (SP) 
dr dr 
and 
do N2 
” oe 
R (Fyr ile (2 


(B-16) 


z SES | =D: 
; n2 5 ERESe Å 
al +r? (5) + sine (7) sin?9 
Å dr dr 


C and D being integration constants. From (B-15) and (B-16) we get 


ye EGEDE 8 (ey BH! DE RER ENDER É 
| Hr (€3 + sin? He i ) Yale i (B-17) 


SS 


| JEG) || r K (7) 


On a ray emerging from the point P we consider a variable point P'. The coordinates 
of P'" with respect to a polar coordinate system with the axis passing through P are 
denoted by (rr, AX, Ap). In the spherically symmetric case, the rays must obviously 
be plane curves, and therefore Ap is constant along any ray through P. Using the 
coordinates (rp7”, AW, Amp) instead of (r, &%, p) in (B-15), (B-16), and (B-17), we there- 
fore find that C = 0 and 


dåAÅv D D (B-18) 


festen rr ER D DD 
drp FALK (Te) /1-( ) P PB 


rpK (Tp) 


and, hence, 
9% 
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ke HE 249 


Der (B-19) 


DE SE NEL 
FK (7) pe dr TEK Gr) 


re IK Ey 


The signs on the right-hand sides of (B-18) and (B-19) apply to the case that rp'<rp, 
if D is chosen to be positive. We shall restrict ourselves to this case, but we remark 
that the signs are to be reversed if rp7>rp. By substituting (B-19) into (B-17), we get 


9 2 
Øg 2 1 i 
—- het SE re sin? (dø VILE RENE (B-20) 
dr 2 i dy 
raket: RE 
A hæs K fø) 
and hence, approximately 
ap" 
(kr as= VK ars LP", (B-21) 
WP erp' ÅG 
where 
Arp 
dr 
c= N ey (B-22) 
AE KO) 


From (B-10) and (B-21) it is clear that, in the integration over P" in (B-14), we get 
essential contributions only when 
ao 


1% 
PAS 


This justifies the expansions which have been made in (B-17), (B-18), and (B-19), 


for in typical cases 
2 


DEER Ad ER ; 2: 
KO % (EMNET DCN 
SK ( SFS ETERN 
GIV 


if re < or < rp and P" lies on the nuclear surface and P lies somewhere in the middle 
of the potential barrier. 

Before we can proceed further we must calculate div es in our spherically sym- 
metric potential barrier. To this purpose, we consider a (somewhat cone-like) surface 
generated by a narrow bundle of rays emerging from the point P, a variable point 
on any one of these rays being denoted by P'. The area which this surface (tube of 
rays) cuts out on the surface through P", which cuts the rays emerging from P ortho- 
gonally, will be denoted by A. For small angles Ad the surface element A is approx- 
imately orthogonal to rr. Using this fact and the formula 
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drp' ER il 
9 ar 
rer) — 
Å BER HR CE) 


which follows from (B-18) and (B-19), we get 


lA ; | 
ER  —= 2 SEG 555 me É: Fa i å Å . 
ES [1 DP; 
r%, K(rp' | ESSE 
He, Bær (7) 


Furthermore, according to Gauss” integral formula, we have 


dÅ 1 dÅ 


Æw 


SEE ES 


div es 
(ther pontERd) 
and, hence, 


div es Sv | - fan | å (B-23) 
(at the point P'") 2 i Bd La 
[BOT TEK (CP | EET 
BEL RA En r? K(r) 
sf RÅ BRS small" 
From (B-23) we easily find that 
PINE [ P' | 1 
(Å ds) exp i 1 divesds, — RER (B-24) 
we | Pi 2 il 
EP E re re K (r?)N mr ER 
Jr? GG: 


if A9 is small and P; is chosen to lie very close to P. From (B-10), (B-21), and (B-24) 
we then get the approximate formula 


eee TE HV: — (49)2| KA 
ZA D GOSS reed RGD EDNBVE GE AJ (B-25) 


wherele 1sidefined" by" (B-22). 

We now assume the wave function y(P') at any point P' (rp”, 9", p') on the 
spherical surface r2” = Ro to be equal to Yi, m(9", 9"), and calculate the wave function 
y (P) at the point P (rp, &, p) in our spherically symmetric potential barrier. To this 
purpose, we substitute (B-25) into (B-14). If we then replace, firstly, (es —e): dør” by 
2Råsin(49)d(49) d(49), secondly, (497 by 2(1—cos(49)), thirdly, cos (49) by 
z and use the formula 


VRE (OA DEDE DER 2 SEERE RER CARD RALDVE 


m,n 
n 


we get the approximate result 
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70 
Ro 1 /K (Ro) [| cr SR RR 2 
P)= — BEDER (7) dr ENE AD) (B-26) 
p(P) ÅER Al He | m 
where? 
re i 2 7 i 
ane $ dr= 56 $ ii (-1) (B-27) 
AE 


ji denoting a spherical Bessel function. However, this result cannot be used for very 
large I-values, since the formula (B-14), on which the calculations leading to (B-26) 
and (B-27) are based, is valid only if the wave function varies slowly on the nuclear 
surface compared with its rapid variation in the radial direction. 

The coefficient &, is real and positive, since it is seen, e.g., from the power 
series representing the spherical Bessel function 7; (z) that 7” j; (- ix) is real and posi- 
tive if x>0. For the calculation of &;, according to (B-27), without introducing any 
further approximations, one can start from the formulae 


Ez 
| c 
| Ad 
(been ra c 


and use the recursion formula 


Xl 79 (SU ce 


For an approximate evaluation of (B-27) one can use the formula 
Før me 1 z (El) 
I ÆW -—— SR, —————————————————— 
ii) = sl exp) | po Rå dæ| | 
TT Dr RES See SDR NE 
ge PT ED || Ht =In(|/ ED 4 1/14 1CED)] , 


which is obtained by applying a first order WKB-approximation to the differential 
equation for the function æ j;(— ix), and which is valid if x>>1. If the more restrict- 
ive condition (l (1+1))?/x3 << 1 is fulfilled, the formula (B-28) simplifies to 


SKTED)] 


es £ 1 
tj (Fix) SV ——expræ= 
IE 5 PP) Boe 


(B-28) 


(B-29) 


and from this formula it follows that & is approximately equal to exp Jr 1) 
ifte (I (1 DE <<1. The results of some numerical calculations of & in the different 
degrees of approximation now discussed are shown in Fig. 9. 


£ CHRISTY has also applied his approximation method to the spherically symmetric penetration 
problem and found that the alpha wave function corresponding to the angular momentum / contains the 
I-dependent factor &; given by (B-27) (private communication). 
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Fig. 9. Results of some numerical calculations of &; in different degrees of approximation. The points denoted 
by crosses, circles and triangles (which correspond to c = 0.05, c = 0.10, and c = 0.20, respectively) 
have been obtained from the formula (B-27) without introducing any further approximations. The 
curved line has been obtained from (B-27) and the approximate formula (B-28), whereas the 
straight line has been found from (B-27) and the formula (B-29), which is still more approximate 
than (B-28). In the approximation represented by the straight line, &; is given by the well-known 


expression exp 4- 1 cl (1 + 1)) . 


In the usual treatment of the penetration of an alpha particle of given angular 
momentum / through a spherically symmetric Coulomb potential Up (mr), this problem 
is reduced to the problem of finding the appropriate solution G;(E, r) of the dif- 
ferential equation (A-1) in Appendix A. If the alpha wave function on the spherical 
surface rp' = Ry is assumed to be equal to Yz, m(9”, p'), the alpha wave function at 
a point P inside the potential barrier contains an [-dependent factor corresponding to 
&; in (B-26) and being the quotient of G7(E,r2)/G1(E, Ro) and Go(E,rr)/Go(E, Ro). 
The value which is obtained for this quotient when the strict WKB-approximation 
(A-13a) and (A-14a) in Appendix Å is used for the functions G;(E,r) and G0(E, r) 
will be denoted by »,. In the limit of small values of 1 (1+1)/x?, the formulae (A-13a) 
and (A-14a) can be replaced by (A-15a), and then the above mentioned quotient 
reduces to a quantity which will be denoted by f7, and which is given by 

bre (B-30) 
where 
10= 


Nr. 3 


PET AR oe B eN (B-31) 

a > |be krp N, PRG) 
The quantities &;, $1, and y1 have the same physical meaning, but they are obtained 
by different approximation methods. In Tables 14 and 15, the results of some 


numerical comparisons between these three quantities are shown. It is seen that &1, 
Pr, and y1 are approximately equal up to fairly large values of l. 
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TABLE 14. NABEESLS 
Numerical values of 1/1 Numerical values of &1/P1 and Bi/yi 
for different values of c. for krp/x=0.2, krp/x=0.6 and x=50. 
l ig . l ar/Bi Bilyi 
fore 0,05" |for er = 0:10 fore 0:20 

0 1 il i 0 1 1 

1 0.999 0.995 0.977 il 0.999 1.000 
2 0.996 0.986 0.947 PÅ 0.997 0.998 
3 0.993 0.974 0.930 3 0.994 0.996 
4 0.989 0.966 0.944 4 0.990 0.992 
5 0.986 0.964 1.008 5 0.987 0.985 
6 0.983 0.971 LS 2 6 0.984 0.974 
7 0.981 0.994 12427 7. 0.983 0953 
8 0.982 1.036 1.94 8 0.983 0.935 
9 0.987 A 406 9 0.986 0.905 
10 0.995 1:21 10 0.993 0.869 
sal 1.009 1537 ital 1.004 0.823 
2 1.029 1:61 12 1.020 0.770 
13 18057 1.96 13 1.043 0.709 
14 1.095 2.5 14 1.075 0.642 
T5 1145 u5 12616 0.570 


j After this detailed treatment of the spherically symmetric penetration problem, 
we now consider the penetration of an alpha particle through the anisotropic potential 
V (r) outside the surface of a deformed nucleus which we suppose to be fixed in 
space. We write the equation for the nuclear surface in the general form 


r=R(W), (B-32) 


and we assume that the alpha wave function on this surface is a given function 
po (VW, p'). Assuming the purely anisotropic part V (7) — Up (r) of the potential V (r) 
to be small compared with its spherically symmetric part Up (17), we can in first ap- 
proximation use the same rays as for the spherically symmetric potential Ue (1). The 


reason is that the ray between any two given points P and P" in the potential barrier 
re 


has the extremal property (B-5), which implies that the integral Vi ds along the ray 
TP 


N283 - 73 


between these two points is stationary for small variations of the path of integration 
connecting P and P". According to (B-20), we then have 
2 


ds SK d 
væ WA il + — z rr BR 
dr PÅ LA gs 
ig Ko (7) | om 
fyr KoGe) 


where AW is the angle (viewed from the center of the nucleus) between a point P 
in the potential barrier and a point P" on the nuclear surface, and Kop (r) is the function 
which we obtain from K(r) by replacing V(r) by Uo(r). Hence, we get the ap- 
proximate formula 


op GTP 2 
VEN coon 
P UP i, Hos 
Als rK 
if we neglect terms of the order of magnitude (AK/K7?, where AK = K(r)— Ko (1). 
ap' 
This formula for | K ds can be written 
”P 
ap arp Irp (or p' (49) 
VÆ ds = Ko (7) dr AK (7, 8) de ÅK (2, 9) dr SE 5 
P JR. UR, Ro 2) Ce 
Nee r? K(1;9) 


where Ry is a constant conveniently chosen to be some average value of R(9”). The 
angle dis here a function ofr which is geometrically represented by theray passing through 
ie pomisflFandfP' SIntthellastitermiittistsafeltorreplace re byER or and (7) Dy 
Ko (1). Inthe third term, we can approximately replace the variable angle & by the angle 
ms for the point PS for ir p = Ro | is small compared to Ro, and the most important rays 
go essentially radially in the neighbourhood of the nuclear surface. The third term 
turns out to be more important than the second one, and therefore we may evaluate 
the latter in a rather approximate way. Since AK decreases rapidly with in- 
IT p 
creasing r, the essential contributions to the integral VÆ: dr are obtained for small 
e'R, 
values of r, and therefore we approximately replace the variable angle & appearing 
in the integrand AK (r, 9) by the angle 9” for the point P". Substituting the resulting 
ap" 
formula for |x ds into (B-10), and assuming that P lies so far away from the nuclear 
eP 
surface that the variation in H(P, P") as P' changes slightly is due essentially only 
aP' 
to the change of the integral IV: ds in the exponent of (B-10), we get the approximate 
eP 
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formula 


éR (0) Arp 
EGP ALDJEÆT (BÆR 0) ep) K(r,9”) ar=Å AK(r,9”) ar, (B-33) 

”R, Ro 
where Po is the point on the sphere r = Ro which has the same polar angles (9”, p') 
as the point P' on the nuclear surface, and the index 0 on the kernel Ho (P, Po) in- 
dicates that this kernel refers to the case that the potential is spherically symmetric 
and equal to Ug (1). From (B-14) and (B-33) it follows that, in first approximation, 
the alpha wave function at the point P in the anisotropic potential barrier V (r) is 
the same as if the electrostatic potential were equal to Up (1) (= the spherically sym- 
metric part of V (r)) and the alpha wave function y on the sphere r = Ry were equal to 


| aR (V”) érp | 
wo (9”, p") exp | es Volr= | AK (r,9”) dr [ (B-34) 
”Ro IR 

This approximate result shows clearly how the alpha wave function is affected in a 
very simple and natural way by the variation in thickness of the potential barrier due 
to the shape of the nucleus and by the variation in height of the anisotropic potential. 
After having resolved the expression (B-34) in terms of spherical harmonics, one can 
immediately give the alpha wave function in the potential barrier as a superposition 
of expressions such as that on the right-hand side of (B-26), which we have pre- 
viously shown to be essentially equivalent to 


Ry Gr (E; 72) 


Tp GL(E, Ry) HEADER 


Hence, we get the final approximate formula 


SET Ro Gi (E, Tr) 
d, rl, SER Ø, 
Pe: 2) De r G1(E, Ry) RE) 


. (B-35) 


977 (02 7 IR (VW) er 
x | ly (9, p”) exp hk (17,9) dr — M KEE TAn | Yr Q(W, 9) sind" dd" dy” 
0 Ra 


0 R4 


for the alpha wave function in our anisotropic potential barrier. 


CERN Theoretical Study Division, Copenhagen 
and 
The Institute for Theoretical Physics, 
University of Copenhagen, Denmark. 


References. 


1)"KSALDER, Å. Bo0HR, T: Huus, B. R. MOTTELSON, and Å. WINTHER; a review article on 
Coulomb excitation, Revs. Mod. Phys. 28, 432 (1956). 


2) F. AsSARO and I. PERLMAN, Phys. Rev. 88, 828 (1952). 
3) F. AsaRo, S. G. THOMPSON, and I. PERLMAN, Phys. Rev. 92, 694 (1953). 
A)NESNÅASARO;EFYSTEPHENS, JR, ande ISPERLMAN,»Phys.Reyv: 92,1495 (1953). 
5) F. AsaRroO and I. PERLMAN, Revs. Mod. Phys. 26, 456 (1954). 
6) F. Asaro and I. PERLMAN, Phys. Rev. 93, 1423 (1954). 
7) F. AsaRroO and. I. PERLMAN, Phys. Rev. 94, 381 (1954). 
8) F. AsaRo, F. S. STEPHENS, JR., S. G. THOMPSON, and I. PERLMAN, Phys. Rev. 98, 19 (1955). 
HERSPASAROFandk I PEREMANS PhysyRev-99537(1955): 
F 


10) F. Asaro, F.S. STEPHENS, JR., B. G. HARVEY, and I. PERLMAN, Phys. Rev. 100, 137 (1955). 
11) F. Asaro and I. PERLMAN, Phys. Rev. 104, 91 (1956). 
112) L. C. BIEDENHARN, J.M. BLATT, and M. E. Rose, Revs. Mod. Phys. 24, 249 (1952). 
12) J. M. BLATT and V. F. WEIsskoPrr, Theoretical Nuclear Physics. John Wiley & Sons, 
New York 1952. 
ISYEJSSTBrEarrR PhyssRevt 95711218 (1954): 
14) A. BoxHr, Mat. Fys. Medd. Dan. Vid. Selsk. 26, no. 14 (1952). 
I5)PÅAS Bor PT OSFROMAN, "and" BB: R.-MOTTELSON," Mat: Fys. Medd. Dan. VidsSelskr 297 
no. 10.(1955) ' 
16) A. BoHR and B.R. MotrterLsoNn, Mat. Fys. Medd. Dan. Vid. Selsk. 27, no. 16 (1953). 
17) N. Boxr, Mat. Fys. Medd. Dan. Vid. Selsk. 18, no. 8 (1948). 
18) M. BorN, Optik. Springer Verlag, Berlin 1933. Reproduced by Edwards Brothers, Inc., 
Ann Arbor, Michigan, 1943. 
19) G. H. BRIGGS, Revs. Mod. Phys. 26, 1 and 472 (1954). 
19a) J. P. BUTLER, T. A. Eastrwoon, H.G. Jackson, and R. P.ScHuman, Phys. Rev. 103, 965(1956). 
20) E. U. Condon and G. H. SHoRTLEY, The Theory of Atomic Spectra. Cambridge University 
Press 1953; 
PLVERSES Caeisry, Phys- Rev: 98212057(1955). 
PN RWAT AE DBS ETA DSROBERTS and GEWS PARKERE Bull SA ms Phy ss SocÆæSeries ele 
no. 4, R9 (1956). 
PSYETEJEDEVANEYE Phys: Rev9185877(1953). 
24) G. W. FARWELL, J.E. ROBERTS, and A.C. WAHL, Phys. Rev. 94, 363 (1954). 
25) P. R. Fienns, M. H. STUDIER, J. F, MecH, H. DIAMOND, ÅA. M. FRIEDMAN, L. B. MAGNUS- 
SON, and J. R. HUIzENGA,-Phys. Rev. 94, 209 (1954). 
26) A. M. FRIEDMAN, J.E. GINDLER, R.F. BARNES, R.SJOBLOM, and P.R. Fierns, Phys. 
Rev:102,585(1956): 
27) M. FRILLEY, S. ROSENBLUM, M. VALADARES, and G. BOUISSIERES, J. Physique Rad. 15, 45 
(1954) and 16, 378 (1955). 
28) G. GAMowW and C.L. CRITCHFIELD, Theory of Atomic Nucleus and Nuclear Energy 
Sources. Clarendon Press, Oxford 1949. 
A. GHIorso, B.G. HARVEY, S.G. THOMPSON, and G.R. CHOoPPIN, unpublished data 
(private communication from Professor I. PERLMAN). 
30) R.A. Grass, S. G. THOMPSON, and G. T. SEABORG, J. Inorg. Nucl. Chem. 1, 3 (1955). 
13) L. L. Gornin, E. F. TRETYAKOV, and G. I. Novikova, Report of Conference on Peaceful 
Applications of Atomic Energy, Moscow 1955, Physical-Mathematical Section, p. 226. 


29) 


76 NILSS 


31a) L. L. GoLDIN, G. I. Novikova, and E. F. TRETYAKOV, Phys. Rev. 103, 1004 (1956). 

32) H. J. GROENEWOLD, Mat. Fys. Medd. Dan. Vid. Selsk. 30, no. 19 (1956). 

33) D. L. HizL and J. A. WHEELER, Phys. Rev. 89, 1102 (1953). 

34) R. W. Horr, J.L. OLsen, and L.G. Mann, Phys. Rev. 102, 805 (1956). 

35) J. M. HOLLANDER, I. PERLMAN, and G. T. SEABORG, Revs. Mød. Phys. 25, 469 (1953). 
(Quoted in the text as HPS). 

35a)J. M. HOLLANDER, W. G. SMITH, and J.W. MIHELiIcH, Phys. Rev. 102, 740 (1956). 

36) J. P. HUMMEL, F. S. STEPHENS, JR., F. ÅASARO, Å. CHETHAM-STRODE, JR., and I. PERL- 
MAN, Phys. Rev. 98, 22 (1955). 

37) L. B. MaGnusson, M. H. STUDIER, P.R. Fiernns, C.M. STEVENS, J.F. MEcGE, AA, M: 
FRIEDMAN, H. DIamoNnD, and J.R. HuIzENGA, Phys. Rev. 96, 1576 (1954). 

38) L. B. MAGNussoN, D. W. ENGELKEMEIR, M. S. FREEDMAN, F. T. PORTER, and F. WAGNER, 
Jr., Phys. Rev. 100, 1237 (19535). 

39) L. B. MaGnusson, A. M. FRIEDMAN, D. ENGELKEMEIR, P.R. FiIEnnDs, and F. WAGNER, 
JR., Phys. Rev. 102, 1097 (1956). 

39a)J. F. MecH, H. DiamonDn, M. H. Stunier, P.R. Friends, A. HirkscH, C. M. STEVENS, 
R. F. BARNES, D. J. HENDERSON, and J. R. Hu1zEenGa, Phys. Rev. 103, 340 (1956). 

40) F. F. MoMYER, Jr. and E. K. Hyne, Phys. Rev. 101, 136 (1956). 

41) B. R. MorrezrsoN and S. G. NiLsson, Phys. Rev. 99, 1615 (1955). 

42) I. PERLMAN, private communication. 

43) I. PERLMAN and F. ASsARO, Ann. Rev. Nucl. Sci. 4, 157 (1954). 

43a)I. PERLMAN and J. O. RASMUSSEN, UCRL-3424 (a review article on alpha decay to be 
published in vol. 42 of Handbuch der Physik (Springer Verlag, Berlin)). 

44) J. O. RASMUSSEN, JR., UCRL-2431 (1953). 

45) J. O. RASMUSSEN and B. SEGALL, Phys. Rev. 103, 1298 (1956). 

46) L. D. ROBERTS, J. W. T. DaBBs, G. W. PARKER, and R. D. ErxisonNn, Bull. Am. Phys. Soc. 
Series II, 1, no. 4, R8:(1956). 

46a)L. D. ROBERTS, J. W.T. DABBS, and G. W. PARKER, Nuclear Orientation of Np237, U233, 
and U?35 Nuclei, to be published in Physics Division Semiannual Report, September 
1956, ORNL report no. 2204. 

47) M. E. Rose, L. D. RoBEerTs, and J. W. T. DaBBs, Bull. Am. Phys. Soc. Series II, 1, 
no. 4, R10 (1956). 

48) S. ROSENBLUM and H. TYREN, C. R. 239, 1205 (1954). 

49) G. SCHARFF-GOLDHABER, E. DER MATEOSIAN, G. HARBOTTLE, and M. Mc Keown, Phys. 
Rev. 99, 180 (1955). 

50) G. T. SEABORG and J.J. Katz, The Actinide Elements. Mc Graw-Hill Book Company, 
Inc., New York 1954. 

51) F. STEPHENS, J. HUMMEL, F. AsaRro, and I. PERLMAN, Bull. Am. Phys. Soc. 29, no. 8, 
G5 (1954). 

52) F. S. STEPHENS, JR., F. AsaRo, and I. PERLMAN, Phys. Rev. 100, 1543 (1955). 

33) F. S. STEPHENS, F. ASsARO, and I. PERLMAN, Odd Parity Rotational Bands in Even-Even 
Nuclei, in manucript. 

33a)V. STRUTINSKY, On the Theory of the Alpha Decay of Deformed Nucleus; manuscript 
which will be published as "Letter to the Editor” in Russian in the 2nd Nb of ""Jurn. 
Deore perrskis R956] 

54) G. VALLADAS, THESE, Commissariat å 'Energie Atomique, Rapport C.E.A. no. 483 (1955). 


Indleveret til selskabet den 27, august 1956. 
Færdig fra trykkeriet den 26, april 1957. 


i: 
E… Det Kongelige Danske Videnskabernes Selskab 


Matematisk-fysiske Skrifter 
Mat. Fys. Skr. Dan. Vid. Selsk. 
Bind 1 
(uafsluttet/in preparation) 

1. BRODERSEN, SVEND, and LANGSETH, A.: The Infrared Spectra of Belsehel sym- 
Benzene-ds, and Benzene-dg 19563 SE EPE NAR RER AERO EA SUE HR EEREEGERE SR 

2. NorLunn, N.E.: Sur les fonctions hypergéométriques d'ordre. supérieur. 1956 .… 1 
Froman, PER OLor: Alpha Decay of Deformed Nuclei. 1957................... 2 


Sa 


On direct BR to the agent of the Academy: EJNAR MSN Ø s, 


in 8vo. The Meddelelse! and the Skrifter only differ in size, not at all in in, Sabjelte treate 
Papers with large formulae, tables, plates etc. will generally be published in the "Skrifter 
4to. Since it is due to a mere accident, with regard to the subjects treated, if a paper is. pri 
in Meddelelser in 8vo or in Skrifter in 4to, these two series will be sent together automat! ca 
to all subscribers to the mathematical and physical papers of the Academy. i 
For subscribers and others who want to receive only the publications concerning a sing 
group of subjects, there is a possibility of special arrangement with the agent of the Acai em! 
in order to acquire the papers published under one or more of the subjects: Mathem 
Physics, Chemistry, Astronomy, Geology. 
In order to avoid mistakes at registration and quotation the pibliations v will not | 


are the following: 
Vol. 30, No. 1—13, 15, 17, 19-21. — Vol. 31, No. 1-2. 


Æg 


Printed in Denmark 
Bianco Lunos Bogtrykkeri A/S 


